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A numerical procedure for calculating the molecular 
flow parameters from direct physical approach has been 
developedo The procedure has the merits of being accurate 
and flexible In comparison with the Integral equation 
approach^ yet It has the simplicity of Monte Carlo method* 
In particular^ the effects of subliming surface reflection# 
channel wall reflection and molecular collision In 
cylindrical channels are taken Into account in the 
non-obstruction factor calculations* The same procedure Is 
extended to compute the mass and thrust non-obstruction 
factors for conic channel s* The resultant molecular 
non-obstruction factors and the evaporation coefficients are 
employed in the transient thermal analyses of an erecting 
shell and a subliming microthruster proposed for orientation 
and despinning of a small solar probe* Transient numerical 
solutions of zero to two dimension heat conduction equations 
for both the non-linear and moving boundary conditions are 
Involved In the analyseso The geometry of a microthruster 
utilizing the subliming material and radioisotopes is 
optimized by the above numerical procedure* Alpha-emitting 
radioisotopes# in addition to supplying the heat and thrust 
for a mlcrothruster# are found to Improve the performance of 
an erecting shell designed to damp out undesired 
precesslonal motions by significantly increasing the optimum 
spin rate* Finally# the applications of the erecting shell 
and the subliming microthruster to the attitude control of 
the MIT®s Sunblazer appear very attractive In several 
aspects* Their future roles In the attitude control of 
small solar probes are discussed in some detail* 
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CHAPTER 1 


SUMMARY 


The present study was spurred by the fascination and 
the challenge of designing a passive attitude control scheme 
for small solar probes of the Sunblazer typCo The 
performance of the solar cells employed for electric power 
and of the antenna used for communications can be Improved 
If such small probes can be oriented and the desired 
orientation maintalnedo The usual orientation directs the 
solar cells toward the sun and antenna toward the earth 
during superior conjunction^ to obtain maximum power and 
antenna gaino 

The problems of passive orientation of solar probes 
differ from that for passive earth satellites^ since the 
solar gravity gradient and the solar magnetic field are too 
insignificant to be of practical use In attitude controVo 
As a result^ solar radiation becomes the only directional 
energy source available for solar probeso The approximately 
collimated solar radiation will induce an asymmetry of 
temperature distribution around a spinning spacecrafto 

The two principal objectives in this work were (1) 
study the applications of radioisotopes and subliming 
materials to magnify the effect of this temperature 
asymmetry so that the spacecraft can be forced to orient 
toward the sun rapidly and reliably^ and (2) the 
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optimization analysis of a subliming microthruster for use 
In the initial despInnIng the spacecraft after Injection and 
In long term attitude control « Other required studies 
related to these objectives were? C3) the application of a 
modified Monte Carlo technique to Knudsen or molecular flow; 

the various numerical solutions (involving from zero to 
two dimensions]) of the parabolic partial differential 
equation for heat transfer with non"linear and moving 
boundary condltlonso 

.,L L Jhfi rma.l Aoatols.., 9f , .f ,rsc,t.l. .§h.e.,U . 

A spinning cylindrical shell under the solar radiation 
will develop an asymmetric temperature distribution around 
the shell as shown in Figure Id due to the fact that the 
solar Intensity on the shell surface Is not uni form”" rather 
It Is nearly a half-rectified sinusoid at any point on the 
shell o The transient temperature solution can be 
approximated^ depending on the problem^ by Ignoring the 
circumferential heat conduction effect^ or the radial heat 
conduction effect^ or by assuming no thermal conduction 
effect at al 1 o The results obtained Indicated that for a 
thin metallic shelly the zero-dimension^ or Strip 

Approximation^ is as good as the 6- dimension^ or Thin Shell 
Approx I mat ion o However^ If thickness Increases^ the 

smoothening effect due to the circumferential heat transfer 
also increaseso For such cases^ the Thin Shell 

Approximation is recommendedo For a non-conduct I ve shell of 
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any significant thickness^ the Thick Shelly or r"difnensIon<, 
Approximation is needed to account for the significant 
temperature gradient across the shell thickness » For the 
types of design Investigated^ the Strip Approximation is 
suitable for a bare shell and radioisotope-heated shelV.? but 
the Thick Shell Approx imat ion is definitely required for a 
subliming material coated shell o In addition to the above 
transient numerical methods^ first order analytic 
approximations to the quasi -steady state solutions were 
Obtainedo These analytic solutions have an inherent 
uncertainty of about twenty-five percent^ but they are very 
useful for parametric analysis^ In particular for the 
determination of the optimum spin rate at which the 
resultant force normal to the solar radiation Is a maximumo 
The steady state temperature distribution typically 
resembles a skew sinusoid whose maximum lags behind the 
maximum solar intensity as shown in Figure lolo A net 
thermal reradiative force results^ which has a component 
normal to the solar directiono This normal force component 
can be employed to erect the spacecraft^ loeo<^ reduce the 
precession cone angleo The optimum spin rate Is obtained by 
maximizing the following Integral for the normal force with 
respect to the spin rate 





cos€d6 


(lol) 


where T(c«)^®) is the steady state temperature distribution 
for spin rate CO as shown computed for Figure lol by the 
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numerical method or by the following approximate analytic 
expression^ 

[l+M(w)sln(©+f(w))] (lo2) 

The temperature distribution for any spin rate can be 
characterized by three parameters^ the average temperature 
TCc^)ov 3 ^ the amplitude factor^ and the phase lag<(, 

^(w}e Both the analytic and the nunerlcal solutions to the 
above optimization have been obtalnedo The results ares Cl) 
for a bare thin shelly the optimum spin is about a few 
tenths of a rpm (which Is so slow as to be marginal from 
dynamics point of view)j (2) for a radioisotope-heated 
shelly the optimum spin Is about a few rpm (which Is quite 
feasible); (3) for a subliming coated shelly the optimum 
spin is only a few hundredths of a rpm (which is too low to 
be operative)o The corresponding maximum normal forces per 
unit height of the shell are Cl) 10 dyne/cm for a bare 
shell; (2) about 20| higher for a heated shell; (3) in the 
order of few hundred dyne/cm (note that even at 3 rpm^ which 
Is well above the optimum spin^ the normal force for a 
subliming shell Is about 1 dyne/cm) « 

The transient behavior of these devices have been 
analyzed in detail » The transient temperature maximum^ 
minimum and phase lag as well as the normal and parallel 
farces generated all exhibited damped oscillations such as 
shown in Figures lo2^ lo3# lo^<i. lo5 and lo6o The transient 
usually lasted no more than an hour provided the attitude 
was fixedp Hot^ver^ during actual stabilization^ the 
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attitude Is changing^ therefore^ the temperature 
dtstrlbutfon Is expected to be In transient constantly o 

Studies on Subi lining Mlcroth ruste r 

To optimize the perforniance of an erecting shell as a 
passive control system for a solar probe It Is necessary to 
bring the spacecraft to about the optimum spin rateo 
However^ the original fast spinning spacecraft cannot 
reliably be despun to the desired spin range by a Yo-Yo^ 
which Is a simple device designed to absorb the spin angular 
moBientum fropn the spacecraft « 

As a res«lt<> a subi Iroing ralcrothruster Is proposed for 
further despln« The general configuration Is shown in 
Figure lo7o The sublimation heat Is supplied from 
radioisotopic heat generated by a thin coating of long'” lived 
radioisotope at the bottom of the cup and from a small 
concentration of shorter-lived radioisotope mixed In with 
the subi imlng material « 

The temperature distribution is very Important so that 
the transient sol ut Ion of a heat conduction equation Is 
necessary^ The simple cylindrical cup design was analyzed 
by both explicit and the Implicit difference schemes to 
obtain the transient temperature dt strlbutiono The results 
agreed very well o However^ because of the moving boundary 
condition (due to surface subi I mat Ion) ^ the implicit method 
appeared^ super lore For a sophisticated design such as the 
chopped cone microthruster the detaned te^parature 


18 



distribution required solution of the two~dlmenslon 
transient heat flow equatlono An ADI 
CAlternate''DlrectIon-‘tmpl IcI t) technique was employedo Due 
to the non~ rectangular geometry^ there was a nutnerlcal 
Instability limitation^ which however was not too serlouso 
In general^ for a rather refined mesh net^ a time step of 
one seGond Is quite stableo 

Emphasis during the microthruster study was placed in 
the optimization of performanceo The total thrust comes 
from the summation of normal raoraentimi <rav*cos8) of a11 the 
exit molecules per unit timex. By varying the nozzle cone 
angle and Its lengthy the (cos6) factor can be approached 
toward unityo The numerical technique developed for the 
molecular flow study and mentioned In Section 1«8 was 
adapted to analyze this problem. The results are shown in 
Figures l.S and 1.9. 

To Improve performance further/ the (m*v) factor offers 
attractive possibilities. If the nozzle wall Is heated 
radlolsotopically/ the average speed of the molecules and 
hence the thrust can be Increased because the molecules 
adsorbed by the wall will receive the thermal energy from 
the heated wall. This fact Is shown In Figure 1.10. In 
addition/ performance could In general Improved by molecular 
dissociation/ but this phenomenon was not taken Into account 
In this study because at the temperatures considered here 
the dissociation effect Is negligible. 
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A subi Irnat ion ratcrothruster can produce thrusts of a 
few hundred dynes^ but Its effective life Is qui t^ 1imlted» 
By providing a thin layer of alpha-emt ttlng radioisotope 
coating to the bottom of the thruster cup^ a long duration 
thrust^ of the order of tenth of a dyne^ can be provided 
over the life of the radIoIsotope« 

Problems of mol ecular flow ere Implicitly related to 
the present study through two dimensionless quantities (1> 
the evaporation coefficient^ which Is the ratio of molecules 
actual ly subl Imed from a surface to the theoretical value/ 
and (2) the non-obstruction factpr/ which Is the ratio of 
rarefied gas (mass or momentum) flow rata through a walled 
channel to the rate for the case no channel walls are 
presents The classical approach of solving Clauslng's 
integral equation ( 16^ 17) has limited solutions to simple 
channel geometries and idealized molecular behaviorso 

In principle^ a straightforward Monte Carlo method 
could be applied to this problem. However^ in traoing the 
history of each molecule/ three or four random numbers are 
requi red tp. detarmlne its orientation as wel 1 as Its energy^ 
and depending on both the channel geometry and the assumed 
possible Interactions/ many more random numbers are needed, 
if an adequate number of case histories Is to be considered/ 
good accuracy I s not obtainable wl thin practical computation 
1 1 mes « The comp 1 ex I ty I ncr eases rap I dl y wi th I nc reas I n g 



channel 1ength« This difficulty Is very similar but not as 
complicated as calculation of the deep penetration of the 
nuclear radiation through a thick shield. 

To deal with the problem^ the fpl lowing modified Monte 
Carlo method was developed. Molecules entering a channel 
are assumed to flow from a finite number (a few hundred} of 
point sources from which the molecules are emitted with a 
Maxwel l tan di stribution. The total angular di stribution 

from each point source Is divided Into many groups each of 
which is characterized by a polar ahgle^ an azimuthal angle^ 
and an ampl I tude according to the distribution function 
(note expected values are used). The fate of every 
molecular beam bundle Is followed and Is recorded. Those 
molecules adsorbed by the channel wall require further 
tracing because at steady state they will be reemItted as 
fast as they are adsorbed, A spherical cosine is assumed 
for these reemItted molecules. Their fate Is found by 
considering the probabilities of the various fate of each 
molecule reemItted from any position In the channel. 

The mathematical variance of this method Is not 
analyzed here^ but In comparison with results obtained by 
other methods for two simplified assumptions (JL£/,LZ.)^ this 
method gives excel lent agreement; see Tables 1,1 and 1,2, 
The main advantage of the new method 1 les In Its flexibility 
both to account for complex raolecular Interactions and 
reflections from surfaces and to be adapted to channels of 
non- uni form cross section. 
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The results for the subltinatlon surface reflection 
effect^ the molecular collision effect^ and the conic 
channel effect are shown in Figures lolly lol2 and lol3o 
Using these accurate channel obstruction effects^ the raw 
data from mass*°flow experiments (JJ,) can be Interpreted to 
produce accurate values for the evaporation coefflcients» 
The result of a sample calculation Is given In Figure lol4o 
The evaporation coefficient is essential In determining the 
heat and mass transfer on a system using a subliming 
material such as the erecting shell and the mlcrothrustero 

1.5 AdpI 1 cat Ions to Sunblazer 

As a dynamical model a small solar probe such as 
Sunblazer is analogous to a spinning top under gravity with 
the exception that the solar pressure acts like a negative 
gravity torque. Assuming the nutation or the wobble has 
been taken care by a viscous damper^ the attitude control is 
reduced to the problems of despin and precession damping. 

Two pal rs of mi crpthruster cups loaded with about one 
pound of subllmlhg material could despin the spacecraft from 
200 rpm to few rpm In less than a day, Howevery If the 
ralcrothruster is used to supplement the action of a Yo^-Yo, 
the requi rement for subliming material is reduced to a 
fraction of a poundo The long-term thrust from a coating of 
alpha emitter at the bottom of the cup is about a thousand 
times smaller than the thrust by subllmatlono Unless the 
roicrothruster Is mechanically turnabley the residual torque 
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of about 1 dyne-cra will continue to change the spin rates 

As the rate of spin Is reduced to the optimum region^ 
the erecting torque generated by an erecting cylindrical 
shell reaches a peaks For a shell coated with subliming 
material; even at the non-optimum spin of a few rpm^ the 
resultant erecting torque Is enough to erect the spacecraft 
In a few hours s For a bare shell; at the optimum spin rate 
(0s37 rpm); the erecting period Is In the neighborhood of 
two weeks s Finally for the radl.olsotope-heated shell; the 
erecting torque Is only slightly larger than that of a bare 
shell; but the optimum spin Increases several folds Since 
the rate of erection Is approximately proportional to the 
ratio of erecting torque to spin rate; the erecting period 
Increases to about five weekSs 

Other attitude control procedures; using an erecting 
shell and/or mlcrothrusters; can be designed for small space 
probesc They may be completely passive; semi -pass! ve; or 
even active such as by synchronlzlns the turning of the 
micrpthrusters as the vehicle rotates so that a desirable 
averaged precession torque and erecting torque resul ts 

(1> The asymmetric temperature distribution around a 
cylindrical shell Is very suitable for generation of 
erecting torque to stabllze solar probess Since the 
transient period Is short compared with the despin period 
following exposure of the shell to the solar radiation; Its 
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effects were Ignored for the present studVo In detailed 
dynamics analysts^ the transient effect always exists during 
the stabilization period because the attitude of the 
spacecraft is continuously changingo 

(2) A subliming material coated shell is by far the 
fastest erecting device but is a1 so the most short*^! i vedo 
The radioisotope-heated shell is capable of the least 
erecting rate but It Is long-lived and can be operated at a 
higher spin rate where disturbing torques are not likely to 
dislodge Its orientationo 

C3) The numerical technique developed for the 
non-obstruction factor calculations of molecular flow is 
very flexible^ accurate^ and can include complex physical 
phenomena o It should also be applicable to the problem of 
calculating neutron and ganvna ray streaming through gaps or 
beam-ports o 

(4) The modified Monte Carlo technique seemed to be 
very appropriate for analyzing the vapor flow for the 
microthrustero The technique eliminates many unnecessary 
complex mathematicso The result clearly pointed the ways to 
improve the design of a microthrustero 

(5) The applications of the above devices^ with proper 
feed-back control mechanism^ can be extended to a wide 
variety of satellites« In particular^ a pair of 
microthrusters can produce torques in all three body axes to 
control the attitude of a 1 ight spacecraf t» 
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FIG 1.1 STEADY STATE TEMPERATURE 
DISTRIBUTION AROUND A BARE SHELL 
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ANGULAR POSITION , (9. deg 
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FIG 1.2i TRANSIENT MAXIMUM 
TEMPERATURE OF BARE SHELL 
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FIG 1.3 TRANSIENT MINIMUM 
TEMPERATURE OF BARE SHELL 



o 

Ln 


o 

o 


O 

tn 

«M 


o 

o 


o 

in 


o 

o 


o 

in 


o 


27 


TIME t.sec 



FIG 1.4 TRANSIENT PHASE 
LAG OF BARE SHELL 
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FIG 1.5 TRANSIENT NORMAL FORGE 
COMPONENT FOR BARE SHELL 
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PARALLEL FORCE, F,, 10"" dyne/cm 


FIG 1.6 TRANSIENT PARALLEL FORCE 
COMPONENT FOR BARE SHELL 
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FIG 1.8 OPTIMIZATION OF CONE ANGLES 
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FIG 1.9 

EXTRAPOLATION OF THROAT AREA 
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L/R RATIO 



THRUST NON-OBSTRUCTION FACTOR, 


FIG 1.10 

HEATED NOZZLE EFFECT ON 
THRUST NON-OBSTRUCTION FACTOR 
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MASS NON-OBSTRUCTION FACTOR, 


FIG 1.11 SUBLIMING SURFACE 

REFLECTION EFFECT 



EVAPORATION COEFFICIENT 





NOfshOBSTRUCTION FACTOR. 


FIG 1.12 MOLECULAR COLLISION EFFECT 
ON NON-OBSTRUCTION FACTOR 
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FIG 1.13 

NON-OBSTRUCTION FACTORS 
vs L/R FOR 45' CONE 
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LENGTH TO RADIUS RATIO, 



FIG 1.14 

EXTRACTION OF EVAPORATION 
COEFF FROM MASS FLOW DATA 
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TABLE lol 


NON-OBSTRUCTIGN FACTOR OF CYLINDRICAL TUBES 
WITH REFLECTION GGEFFICIENT EQUAL TO ZERO 


L/R 

(THIS WORK) 

K 

(CLAUSING) 

Q 

( DEMARCUS) 

0.5 

0.80170 

0.8013 

0.80127 

1.0 

0.67192 

0.6720 

0.67198 

1.5 

0.58165 

0.5810 

0.58148 

2.0 

0.51486 

0.5136 

0.51423 

3.0 

0.41974 

0.4205 

0.42006 

4.0 

0.35619 

0.3589 

0.35658 

5.0 

0.31074 

0.3146 

0.31053 

6.0 

0.27823 

0.2807 

0.27547 
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TABLE 1,2 


NON-OBSTRUCTION FACTOR OF CYLINDRICAL TUBES 
WITH REFLECTING WALLS BUT NON-REFLECTING ENDS 

For L/R«2,0 


r 

Q 

K-M 


(DEMARCUS) 

(THIS WOl 

0,0 

0,51423 

0,51486 

CM 

o 

O 

0,58247 

0,58237 

o 

o 

0,65890 

0,65651 

0,6 

0,74690 

0,74163 

0,8 

0,85412 

0,8 4483 

1,0 

1,00000 

1,00000 

K — 

SEE REF, (li) 


Q — 

SEE REF, (17) 
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CHAPTER 2 


INTRODUCTION 

2.1 InfePQductory Remarks 

The Sunblazer eatel 1 1 te which Is being designed 

and developed at the Center for Space Research of Mo I « To 
will be used to map the charged particle density In the 
ylclnlty of the sun and study the solar corona and magnetic 
fteldo This small spacecraft Is to be powered by solar 
cells which supply an average power of about 20 watts of 
electricity to, a multi -frequeriicy pul sed-transml tter and 
other e1ectronlos<> The transmitter Is designed to radiate 
several pulaed signals simultaneously which wpul<^^ arrive at 
receivers on Earth delayed In time by different amounts due 
to the dispersive nature of the Interplanetary charged 
part i ol as 0 

An iroportant characteristic of this spacecraft is Its 
low weight so tha^ li^cgpenslve Scout rockets can be employed 
to Inject sevafal probes Into hel locentr Ic orbits per yearo 
In do I ns so^ systematic measurements of corona electron 
dens I |y fluctuations could be made « The success of the 
spacecraft depends upon its ability to orient Itself 

continuously to obtai n steady radi atl on to the solar cel 1 S/ 
for, proper thermal control of the spacecraft and for mazlmura 
transmitter antenna gef ho The payload limitation of the 
Scoutrtype rocket precludes the possibility of using 

conventional active attitude control devices such as gas 



jets or reaction wheels which require active on-board 
feed-back control systeoiso However reliability 

cons (derations dictate a system with a minimum of moving 
parts (pass i ve or semi opass i vel while ml ssl on requi rements 
call for a lightweight system which has a relatively fast 
response-time to shorten the period of stabllizattonft The 
purpose of this work is concerned with the design of 
suitable dev ices of this type and the related physi cal 
probl ems« 

By retrograde injection^ ioe,# launching the rocket 


directly 

oppos i te 

to Ear 

th*s orbital velocity# 

the 

spacecraft 

will fall 

into a 

heliocentric orbit with 

i ts 

per Ihel ion 

1 ess than 

1 AU ( i> e 

0 # on an Inferipr orbit) » 

A 

typical 

Sunb later 

orbit In 

an ecliptic projection 

In 

rotating 

and sun-or 1 ented 

coordinates# such that 

the 

Earth’^Sun 

1 i ne stays 

f ixed^ is 

shown in Figure 2» 1* 

The 


information obtained at superior conjunction Is of greatest 
interest^ and consequently the time to super ior conjunction 
(occultatlon of spacecraft by the Sun) should be made as 
short as possibleo Howeyer^ the shorter the superior 
conjunction time or the mission tisie^ the larger the booster 
requirement and the closer the perihelion Is to the sun* 
The small perihelion Introduces probl as>s of both spacecraft 
thermal control and the solar-cell degradation problemso 
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The transnil sslon of pulsed signals at superior 
conjunetlon is parttcularly important for deterinf ning the 
electron densttles In the solar corona. Unfortunately^ a 

communl cat ion blackout region which is a function of 

solar activity as well as the transmission power^ exists 
near occultation, in order to minimize the blackout period^ 
it Is highly desirable to point the spacecraft directly at 
sun so that during superior conjunction the solar cell panel 
and antenna are facing the sun and earthy respectively^ 
which will cause both the electrical power and the antenna 
gain towards the earth to reach a peak. The present study 
Is concerned wi th passive attitude control devices suitable 
for this type of solar probe. The performance of the 
proposed devices are discussed in Chapter 4 and Chapter 5. 

Under direct injection^ the probe will assume a 
superior orbit (with the apehellon greater than 1 AU). 
Severval difficulties arise to affect the attitude control 
directly or indirectly. Owing to the Increasing excess 
velocity requirement for a superior orbit injection^ the 
payload has to be kept small. On the other hand^ because of 
the larger orbit and the longer mission time^ the solar 
Intensity decreases while the transmission power and the 
orientation requirement Increase rapidly with distance. 
Consequently eonceptual designs which depend upon solar 
radiation for power generation and orientation probably can 
not be applied for superior orbit missions. Some of the 
possible solutions are discussed In Chapter 6. 
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2 >5 Spacecraft Dvnamles Consldei^atlons 

2.3 ol General 

For a spln-stabf 1 ized spacecraft/ the ipportant dynamic 
potions consist of hutatIon> precession and spino 
Nutjatlon^ soraetlmes referred to as the force-free 
precession^ Is that motion of a spinning body which deviates 
from uniform spin about a fixed axis in the absence of 
disturbing torciues. Precession is that motion of a spinning 
body which results from an external torqueo 

The present study does not consider nutation damping 
because a simple^ light weight/ yiSGOus liquid damper such 
as an annulus partially filled with mercury and mounted 
concentric with the spin axi s, at a point off the axial 
center of gravity can effectively diss I Pete the nutat lonal 
energy (4)* This Is true only if the spacecraft is 
spinning about the axi s of largest moment of inertia (^) 
because the resulting conversion of kinetic energy into heat 
through viscous damping must eventually bring the body into 
a minimum-energy stateo (Note the terminal rotational 
kinetic energy Is Ta|(»)*/2» where H Is the constant 

angular momentum; for minimum the moment of Inertia I 
must be maximum) o Unfortunately^ this nutation damper is 
Incapable of removing the processional motion because 
viscous damping usual 1 y generates no external torque to 
change tJ>e angular momentum of the spacecrafto 
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Assuming the nutat tonal mot f on can be rapidly damped 
out by an annular liquid damper/ the attitude control 
problem Is reduced to that of spin control and precession 
dampings 

A complete attitude control system must have the 
capability to produce torques about all body axes of the 
spacecrafts These are the spin torque/ the precession (or 
restoring) torque/ and the erecting (or precession damping) 
torques The positive sense of each torque Is defined as 
shown in Figure 2s 2s 

2»3«2 Restoring vane 

Falcovltz (J,) has shown that a sun-orbiting spacecraft 
can be designed to take advantage of solar pressure (loe« 
photon momentum) for at 1 1 tide control by an ingenuous 
arrangement of its external surface — in particular^ by 
corrugating the optical coating on the surfaces of the 
restoring vane or the solar sail for generation of torques 
about all three body axess Unfortunately/ this technique 
can not produce a large non-conservative erecting torque to 
dissipate the kinetic energy stored in unwanted precessional 
mode. Furthermore/ the effectiveness of this technique will 
be decreased by degradation of the surface radiation 
properties of the coated material under long-term intensive 
solar radiation exposure <> The latest Sunblazer design (i.) 
uses four stepping motors to rotate four big triangular 
vanes in order to control the resultant average torques 
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about all three body axes. The torques^ whioh are time 
varying in each spin cycle/ result froin the vane surface 
reflection. The disadvantages of this design are: (1) 

necessity of deploying four large vanes when in orbit; (2) 
no independent control on the magnitude of each torque; (3) 
the capability of attitude control depends on the proper 
functioning of the stepping motors; and (4) out of phase 
electric power requirement/ i,e, more power Is needed to 
drive the motofs when the solar ceils produce jess than 
normal power o Nevertheless/ a restoring vane Is desirabie 
because of Its contribution to the Inherent stability of the 
spacecraft by positioning the center of pressure behind the 
center of mass. The resultant torque will force the 
spacecraft to have e retrograde precession/ i,eo/ to precess 
in a direction opposite the spin direction, 

2 , 3 , 3 E rec 1 1 n g $ he 1 1 

Consideration of gyroscopic action suggests that an 
erecting torque can be generated only if there is a force 
acting on a spinning body normal to both the spin axis and 
the axis of precession, A very clever passive device (2> is 
to simply place a thermally absorptive thin metallic shell 
in front of the spacecraft. The approximately parallel 
solar radiation will always introduce an asymmetry of 
temperature distribution around the cylindrical shell. 
Fortunately/ the resultant rerad i at I on force has a normal 
component to erect the spacecraft or to decrease the 



precesston cone angle as shown In Figure 2«3, in principle^ 
this passive erecting torque should eventual 1y damp out the 
precession and should orient the spacecraft to point at the 
sun dlrectlyo However^ because of the sinali orbital angular 
velocity wblch is always normal to the ecliptic^ It is 
impossible to align the spin axis with the sun axiso Hence^ 
even at equilibrium, the spin axis will point above or below 
the ecliptic plane depending on the spin di recti on » 
Typically, the equi librium canted angle is of the order of a 
few degrees (!)«, 

The serious drawbacks of this simple design ares (1) in 
order to maximize the radiative erecting force, the shell 
has to be very thin and the spacecraft has to be despun to a 
few tenths of a rpm (li.) (As a result of this small spin, 
the spacecraft is very vulnerable to minute disturbing 
forces resulted from the meteorite Impacts, out-gassing. 
Coulomb and Lorentzian interact ions«}; (2) The useful 
reradi ative energy is only a small fraction of the absorbed 
solar energye Also, thermal reradiatlon is the most 
ineffective method of extracting recoil force from a given 
amount of energy to be radlated» 

in Chapter 4, several conceptual techniques involving 
the use of radio! so tope heating and subliming materials 
which are aimed at Improving these shortcomings are proposed 
and analyzed both analytically and numeri cal ly» The 
analytic approach Is primarily for studying the steady state 
behavior in a closed form. The numerical app roach Is to 


4-7 



Investigate the transient and the steady state behaviors 
aGcurately» 

2.«,3o4 Despln raicro thruster 

The residual spin iBinediately after orbit Injection by 
a Scout rocket is about 200 rpm ii) * which is too high for 
stablTizatlon and perhaps for the proper functioning of 
InatrurnentSo A sifnple mechanical device^ the Yo-Yo, is in 
theory capable of reducing the spin to any desirable value. 
It consists of a small mass, m, on the end of a light cord 
wrapped around the spinning spacecraft. As soon as the mass 
ro is released. It will unwind at a constant rate (^) and 
will absorb the spin angular momentum from the spacecraft. 
After the cord Is completely unwound (In a fraction of a 
second), the mass m Is allowed to fly away leaving the 
spacecraft spinning at ^p which Is less than the Initial 
value a)«. The length of the cord determines the final spin 
and it is given by (J^) 

J - Q l(-^ -+ 1 ) 

where R is the radius of the spin body 

I Is the moment of Inertia In the spin axis 

However, for a given cord, the final spin can be known only 

about as accurately as the initial spin Is known. 

Unfortunately, due to the build-up of nutation during 

vehicle deployment, the Yo'^Yo can be safely designed to 

despin the spacecraft only to about one tenth of Its initial 

spin. Hence, further despln is required. In Chapter 5, use 
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of a subi lotion miGrothruster to acccmpl ish the desired 
residudi despin Is proposed and Its performance analyzedo 
The performance of the microthruster is optimized with 
respect to several pertinent varlableso The transient 
solutions of the temperature distrl button In the subi iming 
material are given for both cylindrical cup and chopped cone 
designs 0 Radiolsotppes are proposed to supply the necessary 
heat to the subi Imat Ion process and to generate low thrust 
by expulsion of alpha particles so that the ralcrothruster Is 
capable of the short-term fast despin as well as the 
long-term fine spin control * 

2 . ft Ra ref I ed Gas 0 vn ami cs Cons I de rations 
2,4 Evaporation coefficient 

The evaporation coefficient is the ratio of the actual 
to theorettcal sublimation rates. Experimental values vary 
significantly between different Investigators (9,10,11) » 
Many theories have been proposed to explain this phenomenon 
(9,12,13, 1ft), Burrows attributed the cause to the molecular 
collision (12) , LIttlewood and Rideal blamed the 
discrepancies on the inaccurate thermal measurements (1ft) , 
Some even went so far to suspect the Irregularity of crystal 
surf ace as the main cause (JLL)« Since the evaporation 
coefficient hes a very Important bearing on the design of 
subi I mat Ion devices^ the theory of evaporation coefficient 
and a numerical calculation technique are discussed In the 
second half of Chapter 
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2 . 4 62» Mol ecul ar f 1 ow 

Subllmatfon usually taHes place only at s uch a 1 ow 
pressures; consequently the yapor tnol ecules have a fnean free 
path comparable or even greater than the characteristic 
channel dimens ion o The impedance to the molecular flow 
t;hrough a Ghannel is very irrpprtant for proper determination 
of the evaporation coefficient from experimental data (15.) 
as well as for the design analysis of a subl iraat ion 
microthruster« This flow impedance effect Is accounted for 
by diraensionless quantities such as the Clausing Factor 
(1£)^ Transmission Probabt I ity (lZ) and th this work/ the 
Non-obstruction Factor » These factors all represent the 
ratio of the outlet to inlet raass flow rates through a 
G han ne 1 « C 1 a u s i ns (1& ) ^ b a sed o n geomet r 1 c con s i de r a 1 1 on / 
predicted results for a very s*™PTe cylindrical case of 
col 1 islonl ess gas in a adsoi*PtIve channel » DeMarGus (12) 
applied several advanced numerical techniques to solve 
Cl aosing's integral equation accuratel y and extended the 
results to include the effect of partial specul ar reflection 
of the channel. Mowever/ the Integral equation approach Is 
quite difficult to apply to a conical channel and to use in 
accounting for complicated proceasas such as molecular 
collision and subl imat I on surface refleGtiqn. A modified 
Monte Carlo techhtque was developed In this study and I s 
reported in Chapters 3 and 5/ to simulate the rarefied gas 
flow through both the cyl indrical and the con i cal channel s. 
Because of the directness of this approach/ the method can 
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include various gas*’g;as and gds*‘surface Interactjons as well 
as tha eff®ct of the angulai’ distribution of the exit 
Riotecules on the net thrust. 


5 / 






FIG 2.2 DIAGRAM OF TORQUES 
ON SPACECRAFT 
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FIG 2.3 RERADIATIVE TORQUES 
ON ERECTING SHELL 
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CHAPTER 3 


SUBLIMATION AND RAREFIED GAS FLOW PARAMETERS 
3>1 Introductlen 

The applicatfon of subi imlng naterlals Is suggested In 
this study both as a means of Improving the erecting 
capabilities of a spinning shell and as the fuel for a 
mrcrothruster system for desp Inning and/or stabllzlng small 
solar probeSo In the mlcrothruster the subliming material 
would vaporize In a cup and flow out of the cup to produce 
thrust. The net amount of material sublimed and Its rate of 
flow from the cup are Influenced (l.e.^ obstructed) by 
adsorption and reflection on the cup walls and on the 
surface of subliming material as well as by molecular 
collision between vaporized molecules. Furthermore the 
results of experiments conducted to determine the Intrinsic 
rates of sublimation of various materials are also affected^ 
to various degrees^ by the same phenomena. 

Consequently methods for predicting the effects of such 
*'obstructl ng" phenomena are developed and discussed In this 
Chapter. These methods will be employed In evaluating 
available experimental data on sublimating rates and In the 
application of these evaluated experimental results to the 
design of the stabilizing devices discussed In this report. 
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3»2 Molecular Non-obstruction Factor (Modified Clausing Factor > 
5o2ol Introduction 

The phenomena which affect the rate of flow of 

vaporized material out of a microthruster cup are: 

adsorption on the walls; specular and diffuse reflection 
from the walls and the surface of the subliming material at 
the bottom of the cup; and collisions between the vaporized 
molecules themselves. Since the quantities of interest are 
the rate of material actually leaving the cup and the 
resultant thrust^, the effects of the phenomena mentioned 
will be evaluated in terms of '*non“obstruction factors*'. 

The collective effect of the above physical phenomena on 
non-^obst ruction factors/ will be referred to as the 

Molecular Non-obstruction Factor/ by which the actual rate 
of vaporization from the surface of the subliming material 
can be multiplied to obtain the rate of flow from the exit 
of the raicrothruster cup/ and as the Thrust Non-obstruction 
Factor (or Thrust Factor) which relates the actual thrust 
produced to the thrust of the molecules as they leave the 
s ubl 1 ml ng surface . 

3.2.2 Vapor-wall collision effects 

In a free molecular flow regime/ the mean free path of 
the gas molecule is so large that intermolecular collision 
can be neglected. The rate of flow through a tube is 
limited by collision of molecules with the wall alone. The 
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manner of attacking this problem of flow was first suggested 
by Knudsen (JL^) whose method was to equate the rate of 
momentum transfer to the tube wall by molecular Impacts to 
the difference In the pressure forces at the endso Due to 
his erroneous assuraptlon of the Maxwellian distribution for 
molecules In the channel In calculating the momentum^ the 
free nx>lecular flow equation Is valid only for long tubes 
(proven mathematically by DeMarcus (lZ))o For short tubes^ 
the end effect becomes dominant. Oushman^s formula (18) 
takes this Into account by adding the flow resistance 
through an orifice to Knudsen's long-tube flow formula (U.)# 
analogous to the electric resistance In series. The 
resultant formula Is simple^ but It Is only approximately 
correct because flow resistance does not behave exactly as 
the electrical resistance. 

Realizing the geometric nature of molecular flow^ 
Clausing CJL&) derived a linear Integral equation to express 
the surface col 1 1 si on density function for a cylindrical 
tube, 

^l/(x)-f, (x)+ ^K(y,x)^l/(y)dy 

where the quantity '^(x)dx Is the total surface collision 
density on the wall between length x and x+dx In unit of 
radius. '^,(x}dx Is the contribution due to particles coming 
directly (experiencing no previous surface collision) frcxn 
the Inlet and y is the dummy variable. The kernel K(y,x) Is 
the probability of particles exchange between two 
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dlfferentfal surface areas on the wall of the tube (simnar 
to the view factor in thermal radiation calculations) » The 
desired normalized particle transport rate J is (12) 

d«Js+ A(x)P(x,L)dx 

Js -(#/2)[x^/(x"^^ii)'''*+(x^ + 4)'^*A-2xJ C33; 

The kernel P(x^L) is the probability that a molecule which 
collides with the wall between x and x-t-dx will leave the 
tube via the exit at x"L without further surface coHisionSo 

For ease of obtaining an analytical expression for 
K(y^x) and P(x,L)^ Clausing assumed perfect surface 
adsorption (l<,ep^ no specular reflection) and infinite mean 
free patho The quantity is the fraction of molecules 
streaming through the flow channel directly without 
suffering any surface collisions at all. Clausing obtained 
accurate tabulated results (!£.)> common 1 y known as Clausing 
factors^ as a function of the length to radius ratio without 
the help of modern computing facilities by using asymptotic 
power series expansionso The Clausing factor is defined as 
the ratio of the rate at which gas leaves the outlet of a 
tube to that at which gas enters the inlets in the absence 
of specular reflection, 

in 1955^ DeMarcus (JJ.) investigated the molecular flow 
problem rigorously by starting mathematical formulation from 
the Boltzman equation. Although he obtained essentially the 
same integral equation, he used different numerical 
approaches (^squeezing*, "scissoring^ and variational 
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technietues) to solve this equatIon« According to DeMarcus^ 
the vartattonal result gives the upper limit for the flow 
conductance o He conv>ared his results with Claus log's for 
the case of no specular reflection and concluded that the 
tabulated Clausing Factors are slightly In error because the 
Clausing Factors are In general somewhat higher than those 
of DeMarGUs^ especl al 1 y for long tubes » 

The main drawback of the integral equation method is 
that the formulation of the probability kernels beccmes 
increasingly difficult for non-uniform cross-section 
channels and/or for molecular flow systems with finite mean 
free paths » The method to be discussed adopts the Monte 
Carlo type of approach where the actual physical processes 
are simulated in detail In the form of algorithm for digital 
computer cal cul at i ooo 

3«2»3 Description of numerical approach (TUBNOF) 

The general description of the method proposed for use 
In this study followso For a channel of arbitrary 
cross-section and lengthy as shown in Figure 3al^ the inlet 
is connected to a flat subliming material surface and the 
outlet is connected to a vacuumo (if the back-pressure is 
not zero^ the problem can be divided Into two parts and 
superposition applied to get the net transport.) The inlet 
gas is assumed to be Maxwellian and the inlet area is 
divided into many small pieces^ each of which Is considered 
as a point source. Knowing the angular probability 
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distribution (spherical cosine) for molecules leaving from a 
point source^ the probability distribution function can be 
divided into differential cones as shown schematically In 
Figure 3o2o The fate of the particles contained In the cone 
(Where do they hit the wall? How many? Are they adsorbed or 
reflected?) Is determined by the space curve of the 
intersection of the cone surface with the channel surface 
and the aye raged surface reflection coeffic lento The 

channel Is divided into ring segments so that number of 
particles impinging on each ring segment can be calculated 
from the fraction of space curve intersected by each ring 
segmento The surface adsorption density due to only the 
undisturbed inlet particles as a function of down-^stream 
distance can be obtained by summing up the contributions of 
all the differential cones from a point source and then all 
the point sources in the inlet cross-'sectional areao At 
steady state/ the rate of adsorption will balance out the 
rate of reemissiono Assuming a known angular distribution 
(spherical cosine at the same temperature as the subliming 
surface)/ the probability for an adsorbed particle at x+dx 
distance downstream to get out of the exit can be computed 
by keeping track of this particle until it either leaves the 
channel or is adsorbed by the channel wall»^ Doing the same 
for all ring segments/ an 'escape' probability table is 
constructed which is equivalent to the Clausing's P(X/L) 
kernel o With the surface adsorption density and the 

'escape' probability table so computed/ the molecular 
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nori-obstruct low faetor (also referrecl as CTausIng factor or 
transinl sslon probability) is obtained by numerical 
Integration over the length of the channels The computer 
program (TUBNOF) listing with brief description Is included 
In Appendices Kol and P»8« The appIlGations of this 
computer code to account for various molecular behaviors In 
the flow channel are given below. 

3.2.4 Cal eul at Ions of wall adsorption effect 

The power of the method developed here lies In Its 
flexibility. The method can be employed In the 

consideration of realistic physical phenomenon such as 
Incident energy and angle dependences of gas-surface 

Interaction cross sections as well as to solve for the flow 
In complex geometric channelsy In addition to a few special 
cases of Interest In the present study. To test the 
validity of the method^ the first step was to check results 
computed by TUBNOF against those of Clausing and DeMarcus 
for the case of cylindrical channel/ perfectly adsorptive 
surface/ and Infinite molecular mean free path. The 
tabulated results are given In Table 5.1. The agreement 
with Clauslng^s and DeMarcus's results Is within 1%. 
Consequently/ the method Is believed to be sound. TUBNOF 
is capable of much higher accuracy than thIs/ but a 
significant increase In coiiputer execution time would be 
necessary. 
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3.2o5 Calculations of channel-wall reflection effect 

The general theory for gas-surface Interaction Is still 
not well definedo However^ drawing analogy from neutron 
albedo data and from some recent gas-surface experiments 
C^}/ the reflection (specular and diffuse) coefficient for 
sublimed molecules should be a function of Incident energy 
and angle as well as surface properties and temperature« 
The effect of finite reflection on molecular flow Is to 
Increase the non-obstruction factor, 1 060 , to reduce the 
flow resistanceo By using a weighted average specular 
reflection coefficient, the computer program TUBNOF can 
allow for the specular reflection effect such that whenever 
a beam-bundle (a group of molecules heading In a direction 
of i+d^ and B-i>dB from a point source) hits the wall, the 
fraction that Is reflected speculai'ly will be traced 
continuously until Its Intensity becomes negl I gible due to 
multiple specular reflections or It leaves the channel « For 
a given channel length, the non-obstruction factor, Krr) , 
Increases almost linearly with specular reflection 
coefficient, r, as shown In Figure 3,3p 

A good linear analytical fit to the TUBNOF curve of 
Figure 3p3 for a slightly reflective channel wall Is 

K^(L/R)-K(L/R)aoO-r) + r (3.4) 

where K(L/R) Is the non-obstruction factor 
for r«0 (loep, K(L/R)®the Clausing Factor) 
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The result agrees very well with that of DeMarcus as 
shown in Table 3»2/ which suggests that Smoluchowski *s 
correction 

K5(L/R)»Ka/R>*(1.0+r)/(1.0-r) ( 

Is not applicable for general use as shown In Figure 3o4<> 
In fact, Sinol uchowski *^s correction factor for specular 
reflection is only approx Iinately correct for very Tong tubes 
and for very small specular reflection coefficients® 

5®2<,6 Calculations of subliming surface reflection effect 

For the case where the Inlet molecular source Is an 
equilibrium gas reservoir/ the molecules heading back toward 
the source from the channel will be absorbed in the 
reservoir, I»e«, the source appears to be black to any 
molecules returning from reemIssion or molecular colllslono 
However, when the source is a solid subliming surface at the 
channel Inlet, the solid subliming surface will reflect a 
fraction <1-(X) of the returning molecules, where is the 
evaporation or condensation coefficient (^) ® For a 
collisionless molecular flow, the fraction returning to the 
surface is l-K® Of this fraction will be reflected back 
Into the channel with a probability of leaving the channel 
through the exit of K® This effect will be compounded and 
if the sun of the resulting Infinite series Is considered, 
the sublimation surface reflection should increase the 
non~obst ruction to 
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K* «K/ ( l-( 1-K) ( l-O?) ) 6) 
where K* Is the approximate non*'Obstruct I on factor Including 
the source reflection effect. Consequently^ this surface 
reflection effect Increases with flow channel length (since 
K decreases with length) and should be taken into account In 
the calculation non>obstructlon factors, 

TUBHOF therefore Includes this effect. The (1-Qf ) 
fraction of the returning molecules is assumed specularly 
reflected from the subliming surface with an angular 
direction and intensity which are computed after each 
reflection. The fate of this molecular beam Is accounted 
properly so that no particles will be tallied more than 
once. The result of a sample calculation is shown in Figure 
5. 5. The results of Equation (3.6) agree well with the 
'exact* computer result. 

Through the use of two experiments/ it may be possible 
to determine If finite specular reflection of vapor 
molecules from the subliming materials leads to non-unity 
evaporation coefficients and to determine the magnitude of 
the evaporation coefficient. The experimental technique 
proposed follows, A large equilibrium reservoir with its 
Internal surface coated with subliming material/ with a tube 
penetrating the reservoir wall/ Is used to measure the mass 
flow rate and the non-obstruction factor Kj , Then the same 
tube Is connected to a subliming surface at the same 
temperature as the reservoir and the flow channel length is 
maintained constant so that a new non-obstruction factor 
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can be measured« By using the ratio K;s/K,^a large portion of 
the experimental uncertainty shoul d cancel . The evaporation 
coefficient can be obtained from the curve in Figure 3o6 
which was calculated using TUBNOF or from the approximate 
expression (derived from Equation (3,6) ) 

a »1+(K, /Ka-D/U-Ki ) C3.7) 

which is almost linear as shown in Figure 3 <>6. 

3«2»7 Caicuiations of molecular collision effect 


In the previous sections^ molecules are assumed to 
interact with the surface only/ since the mean free path is 
muoh larger than the characteristic dimension of the flow 
channel 0 This condition does not exist in the present study 
because the sublimation vapor pressure in the temperature 
range of Interest (in the neighborhood of 280°K) has a mean 
free path comparable with the channel dimension» Burrows 
made an estimate of the Intermolecui ar collision effect 
which surprisingly was very signif icant« Burrows in his 
paper (JL2) tried to attribute the cause of the existence of 
evaporation coefficient to the molecular coilisiono It 
should be noted that Burrows' evaporation coefficient 
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where 


K— tube Clausing factor 
a— tube radius 
1— tube length 

%— equil Ibrlura vapor mean free path 
A— effective mean free path 

Is different from the the conventional evaporation 
coefficient by Just the Clausing factor « In other words^ 

the quafitity In the square bracket which accounts for the 
molecular collision effect Is Burrows' estimate of the 
conventional evaporation coefflclento Burrows' explanation 
appears Inconclusive based on the experimental results of 
Sherwood and Uohannes (^) whose sophisticated experimental 
setup should minimize the molecular collision effects 
Nevertheless^ the effect of molecular collision should not 
be ignored^ 

A numerical scheme was developed to estimate the 
molecular collision effect so that the evaporation 
coefficient can be evaluated from the experimental datao 
The first step Is to calculate the molecular density and 
flux (track length density) as a function of radial and 
axial distances (see Appendix 6 for detailed description). 
Only the contributions from the subliming source and diffuse 
reflection from the wall surfaces are considered. For the 
cases of Interest In this study/ the effects of specular 
reflection and molecular collision on density and flux are 
small and are therefore neglected for the sake of 
slropltclty. Due to the lack of knowledge on the dynamics of 
molecular collision/ the direction of motion after a 
collision Is assumed random. One further assumption Is 
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required for the molecular Interaction cross section which 
is related to the effective diameter of subi Itned molecu1e« 
Many well-known techniques (I8\.2p> 31>. 46) of extracting 
molecular size from viscosity^ critical temperature and 
pressure and diffusion coefficient measurements are 
available in the literature. Additionally^ different 
assumed molecular force fields will also yield different 
expressions for the molecular diameter. Since the viscosity 
data Is most readily available, it is used as a basis for 
molecular cross section calculation. The effect of 
molecular collision on the gas flow Is computed In the 
following manner. 

The collision density as a function of radial and axial 
distances is calculated from the result of density, flux and 
molecular cross sections. The solid angles subtended at 
each differential volume by the inlet and exit determine the 
fractions of the collided molecules in the channel heading 
back and streaming out respectively. However, in the 
process of heading back or streaming out, the probability of 
further collisions is governed by an exponential function 
exp(-S/A), where S Is the distance travelled In a particular 
direction and A is the average mean free path in this 
direction. The concept here is analogous to the neutron 
diffusion calculation In a source free medium. Since a 
fraction of the molecules collided was originally streaming 
toward the exit, the net effect Is the difference between 
the fraction heading outward before and after collisions In 
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each differential volumeo The molecular density at each 
position Is computed^ but the angular distribution Is not 
known; In particular^ the fraction of the molecules heading 
out Is not readily obtained^ This fractidn Is assumed In 
the coding of TUBhOF to Increase linearly with axial 
distance and Is constant In radial direction^ The rest of 
the molecules after collision are adsorbed by the channel 
wall in a ratio proportional to the solid angles subtended 
at the collision point by each of the channel ring segments » 
The calculated result for the molecular coirislon effect are 
shown In Figure 3« 7 <, Evidently^ the molecular collision 

greatly Increases the flow Impedance^ 

Also/ since the molecular collision density is 
proportional to the square of particle density^ at low 
sublimation temperature or at low vapor pressure the 

collision effect should be negllglble« The computed result 
Is shown in Figure 3.8, 

in passing It should be pointed out that Burrows* 
expression^ Equation (3.7> probably overestimates the 

molecular collision effect^ because based on Rossmann and 
Yarwood's experimental data ili) * Burrows' corrected 
evaporation coefficients for mercury In many cases were 
greater than unity C which Is not realistic). 



J.S Evaporation Coeff Id ant Theory 


3«3«1 Introduction 

In principle^ all material will sublime from a solid 
phase to a gaseous phase bypassing 1 Iquefactlon/ If the 
transltton occurs below the triple pointo The theoretical 
derivation of the expression for sublimation rate was 
developed by Hertz^ Knudsen/ and Langmuir (19).. The ratio 
of the experimental to theoretical sublimation rate is 
commonly known as the evaporation coefficient« The 
evaporation coefficient/ which is Involved directly in 
calculating the heat and mass transfer rates for the 
subliming microthruster and erecting shell discussed in this 
study/ is therefore an experimental correction factor for 
the Hertz-Kundsen-Langmui r equation (ASlo 

Since the experimental subliming rate is always less 
than or equal to the theoretical rate (^)/ many scientists 
had attenpted to give a satisfactory expianatlono One 
explanation has been that the theory is In error by assuming 
that all collisions with a surface result in condensation 
(the theory assumes that the rates of sublimation and 
condensation in a system at equilibrium are equal and that 
the rate of condensation can be equated to the rate of 
collision on the surface calculated from the known vapor 
properties) o There is also the probabi 1 ity of experimental 
.errorso LIttlewood and Rl deal (lii) went so far to claim 
that evaporation coefficient should be unity If the thermal 
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measurement is performed properl yo Burrows (12) implied in 
his paper that material really sublimed at the theoretical 
rate^ while in actual experin^nts moiecular collisions 
caused the reduction in the subliming rate. Their 
assert ions can not be supported by the recent experimental 
evidences of Paul and Lyon (15) as well as Sherwood and 
donannes (IJ^} whose experiments minimized the molecular 
collision effect and the temperature measurement error. The 
concept of equating the evaporation coefficient to the 
condensation coefficient has been rejected in the past. 
However^ Paul and Lyon's recoil pressure experiment gave a 
favorable support of this concept because the result (based 
on Mtescher's (2) derivation) agreed well with that of mass 
flow. However^ it appears that Miesher's derivation 
requires further refinement which is given in the following 
section. 

3.5.2 Evaporation dependence on reflection 

The detailed theory of gas surface interaction is still 
not well understood. However^ the interactions of vapor 
molecules and their parent solid generally consist of 
diffuse and specular reflections and condensation. The 
condensation process can be understood as a delayed 
reemIssion because molecules condensed will evaporate 
diffusely at a later time and this process is different from 
the prompt diffuse reflection which accompanies the specular 
reflection. In experiment^ the prompt diffuse reflection 
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and delayed reenilsston are difficult to separatee However^ 
In the following derivation of the matheiaat i cal dependence 
of evaporation coefficients on the reflection (diffuse and 
specular) coeff icients, the reflection coefficients as 
functions of surface temperature^ incident angle and speed 
are assumed known o Also the vapor and solid are assumed to 
be at thermal equi 1 ibrjum so that surface cooling or surface 
heating effects are avoidedo Let the total reflection 

coefficient of the vapor molecules on the subliming solid be 

r(«,v^T)=8(e,v,T)+ s(0/v^T) CS'i) 

where &(0^v^T) Is the diffuse reflection 
sO^v^T) Is the specular reflection 
8 is the incident angle 
V Is the incident speed 
T is the surface temperature 

The fraction of the molecules which do not reflect are 
assumed condensed and is given by l~r(8^v^T)o Since at 
kinetic equilibrium^ the fraction condensed would be exactly 
equal to the fraction evaporated^ the total number of the 
molecules evaporated from a unit surface exposed to vacuum 
is given by the Integral 

G = m Ccie r 

-i=0 ■'v=o 

where f (8,v)dvd8“l/2*vdn^slnecos6d8 Is the number 
of molecules per unit time leaving a unit surface 
with speed v and making angle with the normal 
between 8 and d8<> For Maxwellian distribution^ 
the number of molecules per unit volume with speed 
between v and v*dv^ dn^ Is 

C3.ll) 
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The theoretical mass evaporation rate 
HertZ'Knudsen-Langmui r equation is 

✓ 'V/z ^00 _ 

Gm~^ c/8 I dVrCB.'(J) 

Jo Jo 

= Rr/ca'^RTy^^ 

- 0.0 583 Pxr C /W / T ) 

where n is molecular density 
m is the molecular mass 
or is average speed 
R- is vapor pressure in ram of Hg 
M is molecular weight 
T is vapor temperature in°K 
R is gas constant per unit mass 

The evaporation coefficient is then 


or 


a - 


(B. V) dir 


“ I - rce.uT) 

In special case/ when diffuse reflection Is zero. 


evaporation coefficient Is just X-s(6/V/T) 


3«3«3 Ratio of recoil pressure to vapor pressure 


At thermodynamical equilibrium/ the vapor pressure 
the total force acting on the subliming surface) due to 
gas<-surface interactions/ consists of four components » 
fV “Pe +Ps 

where subscript v“-— vapor 

e— — evaporat i on 
d— -diffuse reflection 
s— — specul ar ref 1 ec tl on 
c----condensat i on 

At equilibrium/ the mol ecular condensation rate 
equal to the evaporation rate/ so Pe is equal to Pc 


the 

C3JZ) 


C3.I3) 

the 


(or 

the 

C3l4-j 


is 
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However^ In vacuinn^ no particles will return to the 
subllmlns surface (assutnins roolecular Interaction Is 
negligible); therefore the force acting on the surface Is 
due to Pe alone o Under this circumstance^ Pe Is the recoil 
pressure which Is a physical property of the subliming 
surface temperature and Is Independent of the equilibrium 
condition or the ambient pressureo The expression of Pe In 


terms of Pv Is derived In Appendix E. 


Pe=Rr' 



See, v:t) - rce t) ® + SC&. t) v-cos & 


The barred quantities are the Maxwellian distribution 
weighted averages and 6^ Is the maximum mass flux as defined 
In Equation (3,12), If there Is no diffuse reflection, the 
exp res s I on becomes , 

( 3 . 16 ) 

SCeMT) e- 

Furtherraore, if the specular reflection Is Independent of 
Incident angle and energy, the expression reduces to (using 
Equation (3,13)) 

Pe/Rr = 2 ^ / -SCT)) = 0(/2 ' 

This Is the expression on which Paul and Lyon (JJ.) based in 
their recoil pressure experlmento The validity of their 
result Is dependent upon the condition that the diffuse 
reflection Is negligible and the specular reflection is 
Independent of angle and energy. Consequently, the Pg /Rr 
ratlo from Equation (3,15) Is a function of surface 
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temperature^ but^ In general^ is not a simple function of 
the evaporation coefficient as given by Equation (3«17)o 

3^ 4 Coefficient Calculation Techniques 

3<,4ol introduction 

Precise calculation of the evaporation coefficient from 
experimental data is a very difficult tasko Experimental I y«> 
accurate subliming surface temperature measurement has 
plagued researchers for years « Pau1<. Lyon (JL§.) and others 
(.£) had developed ingeneous surface temperature measurement 
techniqueso However/ in computing evaporation coefficient 
from mas s**f low experimental data/ they generally neglected 
the gas’-gas and gas-surface interaction effects on the mass 
flow rate/ which are quite significant as indicated in 
Section 3«2« In the following sections/ two possible 
techniques for accurate calculation of evaporation 
coefficients from experimental data are discussed. 

3.4.2 Extraction from mass-flow data 

in mass-flow experiments such as described by Paul and 
Lyon (U^)/ both the surface temperature and mass flow 
channel length are a function of time. These two time 
variables as well as the gas-surface Interaction effect have 
to be taken Into account properly before precise evaporation 
coefficients can be obtained. The experimental quantities 
available are the surface temperature and surface position 
as a function of time. The mass flow rate G(t) is 
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(B - 18 ) 


G-(tr) = P 

and the evaporation coefficient 0((t) Is 


act; 


I _ 

G„ct) Kct; 


p 


oix/dil 

GrnCt)!^) 


(3- 19 ’) 


Where Is the density 

Is the rate of surface regression 
GmCt) Is the maxi mum mass flow rate 
K(t) Is the non~obst ruction factor 

The evaporation coefficient as a function of time and 
hence a function of temperature was calculatedo Using Paul 
and Lyon®s experimental data (apparent evaporation 
coefficient) for naphthalene^ the results are shown In 
Figure 3o9 for four cases In which the experimental data 
(bottom curve) were corrected fors Cl) wall adsorption 
effect only or exact Clausing factor correction; (2) wall 
adsorption and subliming surface reflection effect; (3) 
correction (2) plus finite wall reflection effect (r=0ol); 
and (4) correction (3) plus the molecular collision effect o 
Since the higher the non"obstruGtion factor^ the lower the 
corrected evaporation coefficient, it is quite evident that 
Clausing factor correction (1) Is not adequate to describe 
molecular flow accuratelyo The effects of the subliming 
surface and channel wall reflections increase the 
non“obstruct ion factor and hence decrease the corrected 
evaporation coefficient, as shown In Curves 2 and 3„ The 
molecular collision effect as indicated by Curve 4 Is 
responsible for further Increase of the corrected 
evaporation coefficlento Even so the result is lower than 
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the experimental values of Sherwood and Johannes which 
Indicates there may be a residual underestimation of the 
molecular collision effect In the present study, 

3,4,3 Extraction from recol 1-pressure data 

One of the objectives of Paul and Lyon Cl^) was to take 
advantage of the simple relation of Equation (3.17) to 
calculate evaporation coefficients. 

In Section 3,3,3^ the ratio of recoil pressure to vapor 
pressure was shown to be a complicated function of 
temperature as given in Equation (3,15), This ratio# in 
general# Is not equal to the half of the evaporation 
coefficients as given by Equation (3,17), Consequently# the 
recoil pressure method described by Paul and Lyon could give 
erroneous results. Their recoil pressure results did show 
some discrepancies. It is tentatively concluded that the 
recoil pressure experiment is not suitable for extraction of 
evaporation coefficients for the present study, 
Nevertheless# their experimental technique is very useful 
for determination of the thrust non-obstruction factor which 
resembles the mass non-obstruction factor in many ways. The 
thrust non-obstruction factor wi 1 1 be discussed In some 
detail In Chapter 5, 
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no 3.3 REFLECTION COEFFICIENT 
EFFECT ON NON-OBSTRUCTION FACTOR 



REFLECTION COEFFICIENT 

r 
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FIG 3.4 COMPARISON WITH SMOLUCHOWSKI'S 
RESULT ON SPECULAR REFLECTION EFFECT 






MASS NON-OBSTRUCTION FACTOR. 


FIG 3.5 SUBLIMING SURFACE 

REFLECTION EFFECT 



EVAPORATION COEFFICIENT 



FIG 3.6 CURVE FOR INTERPOLATION 
OF EVAPORATION COEFFICIENT 
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EVAPORATION COEFFICIENT, 



NON-OBSTRUCTION FACTOR 


FIG 3.7 MOLECULAR COLLISION EFFECT 
ON NON-OBSTRUCTION FACTOR 





FIG 3.8 SUBLIMING TEMPERATURE EFFECT 
ON MOLECULAR COLLISION 
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FIG 3 9 

EXTRACTION OF EVAPORATION 
COEFF FROM MASS FLOW DATA 



SURFACE TEMPERATURE^ 





TABLE 3pl 


NON-OBSTRUCTION FACTOR OF CYLINDRICAL TUBES 
WITH REFLECTION COEFFICIENT EQUAL TO ZERO 


L/R 

0p5 

K-M 

(THIS WORK) 
0.80170 

K 

(CLAUSING) 

0.8013 

Q 

(DEMARCUS) 

0.80127 

IpO 

0.67192 

0,6720 

0.67198 

1.5 

0.58165 

0.5810 

0.58148 

2.0 

0.51486 

0.5156 

0.51423 

3p0 

0.41974 

0.4205 

0.42006 

4.0 

0.35619 

0.3589 

0.35658 

5.0 

0.31074 

0.3146 

0.31053 

6.0 

0.27823 

0.2807 

0.27547 


TABLE 3o2 

NON-OBSTRUCTION FACTOR OF CYLINDRICAL TUBES 
WITH REFLECTING WALLS BUT NON-REFLECTING ENDS 

FOR L/R-2.0 

r Q K-M 

CDEMARCUS) (THIS WORK) 

OoO O o 51423 O p 51486 

0p2 0 p 58247 0p58237 

Op 4 Op 65890 Op 65 651 

0p 6 Op 74690 Op 74163 

0p8 0p85412 0p84483 

IpO Ip 00000 Ip 00000 

K — SEE REFp (16) 

Q — SEE REF. (17) 
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CHAPTER 4 


TRANSIENT THERMAL ANALYSIS OF ERECTING SHELL 


4.1 Introduction 


The surface temperature of a cylfndrlcal shell such as 
an erecting shell attached in front of a spinning solar 
probe will In general have a non-uniform temperature 
distribution around the shell » It is this asymmetric 
tempera ture distribution which gives rise to a net radiative 
and (or) subliming force to dissipate the precession kinetic 
energy « To analyse the temperature distribution as a 
function of time^ a time dependent two dimension heat flow 
equation of the form (assuming constant material properties) 


5 r' 


U. ^ J- ilL P ii 


r 9r r ^ e*- K dt 


(^> 1 ) 


T(o o, t ) = T(r, 2 Tt , t ) ; 


±1 I 

^ r I r-R+^r 


0 


5T 


5T 1 

* 

aT 1 

r^e 

©so 

ra© 

> 

e - 27T 

5 r 


has to be solved numerically because the boundary conditions 
resulted from the thermal radiation and sublimation are^ in 
general^ non-1 Inear in To The schematic diagram of an 
erecting shell is shown In Figure 4,1, In order to avoid 
numerical instability^ an implicit difference formulation 
has to be used. Consequently the problem Is reduced to 
solving a large matrix equation which is very time 
consuming. In the following sections various assumptions 
are made to simplify the partial differential equation^ 
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under the condition that the results will not be seriously 
compromised. To facilitate comparisons^ the calculations in 
the following sections are based on the Input conditions of 
steady precession (no erecting torque) and uniform Initial 
teirperature. The order of development Isi (1) completely 
neglected the heat conduction effects (Strip Approximation); 
(2) allowed for only the circumferential heat conduction 
effect (Thin Shell Approximation); C3) allowed for only the 
radial heat conduction effect (Thick Shell Approximation), 

,S:tr Ip ^pp rox Imat Ion Me thod 
4 0-2,1 Description of strip model 


The crudest approximation is that the cylindrical shell 
Is so thin that the radial temperature gradient through the 
shell is negligible^ and the shell is assumed to consist of 
strips with Infinitesimal Insulation gaps between strips so 
the circumferential (or tangential) heat conduction is also 
Ignored, The above partial differential equation Is 
slmpl t f led to a set of ordinary differential equations of 
the form 


P 


c 


±I 

^ t 



LF(6i,tKoif-6-<7-T/‘ 


(4.E) 


T^ and 6.,- are the temperature and angular position for 
each strip. Therefore if the shell Is divided Into J 
strips^ there will be J equations to be solved per time 
step. The f(^^t) In the Equation (4,2) is the 
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half-rectif led sinusoid characterizing the solar Input to 
the shell as a function of time and angular position of the 
strip* 

4o2*2 Bare thin shell system 

To maximize the thermal rerad I at I ve effect^ ioe*^ the 
erecting torque^ one of the criteria is to make the shell as 
thin as possible so the solar Input can give rise a maximum 
temperature changeo In this thin shell design^ the strip 
approximation should hold because the radial temperature 
gradient and the circumferential heat transfer would be 
small due to the thin shell thickness* A computer program 
(STRIP) Is written to solve the entire set of Equation 
(4*2)* The description and listing of STRIP are included In 
Appendices H*1 and P*!* 

The angular temperature distribution gradually changes 
from Its Initial uniform shape to skew sinusoid which 
travels around the shell at the angular rate of the spin* 
The Instantaneous maximum and minimum temperatures around 
the shell are shown in Figures 4*2 and 4o3 and the 
instantaneous maximum temperature difference^ a T^^ ^ which is 
the difference of the maximum and minimum temperatures is 
shown In Figure 4*4* immediately after exposure to the 
solar radiation the amplitude of the temperature exhibits a 
damped oscillation character which reaches a steady state 
amplitude after a few spin cycles* 
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The Initial temperature distributions all have the 
damped escilVatlon character which can be explained 
graphically by comparison of the steady state distribution 
with the initial dlstributionp A detailed explanation is 
given in Appendix C where three types of damped oscillation 
are analysed for three classes of initial conditiono These 


are Cl) the initial 

temperature T© 

is 

greater 

than 

the 

steady state maximum 

temper at ure; (2) 

To 

Is 

less 

than 

the 

steady state minimum 

temperature; (3) 

To 

Is 

i n between 

the 


steady state maximum and mini mint temperatures « 

Another Important quantity^ which is plotted In Figure 
is the angular phase lag of the Instantaneous maximum 
temperature from the maximum solar intensity directlono As 
will become clear later^ this phase lag Is a major factor 
In determining the normal component of the reradiative 
force^ and equivalently the erecting torque. 

For fast spin such as the residual spin after the Yo-Yo 
despin operation (see Section 2,3,4)^ the transient 
behaviors will be noticeably different. The most 
outstanding feature Is that the damped oscillation will 
persist for many spin cycles^ but the total transient time 
span Is Increased only moderately. Also^ if the Initial 
temperature Is too high above the equIlibrItBn maximum^ the 
damped oscillation behavior does not exist at all because 
the shell is always losing energy during the entire spin 
cycle. As a result^ the damped oscl 11 at Ion will degenerate 
to a slanted stair shape which continuously decreases toward 
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an asymptotic value, 

4o2,3 Optimum spin rate 

ideally^ to realize the largest norma) force conponent^ 
the phase lag, should be as close to -1^/2 as possible 
and the maximum temperature difference, ATm should be as 
large as possibleo Unfortunately for a given configuration 
at zero spin, is maximum but is zero; and for fast 
spin, AX/i approaches zero but 'vf' Is close to -'Ttllo 
Consequently, there must be an optimum spin rate at which 
the resultant normal force is maximum. Trial-and-error 
solution would certainly give an answer, but it could be 
very tirae-consuming, especially without a prior knowledge of 
the optimum spin region. 

For a re radiative erecting shell, an analytic 
expression for the optimum spin can be obtained by 
maximizing the following integral 

where T(c*>,9) is the temperature distribution 
shown in Figure 4.1. 

The normal and parallel forces resulting from the 
transient temperature distribution of a bare thin shell of 
unit height are presented in Figures 4,6 and 4.7, The 
oscillations have damped out after about 300 seconds for the 
conditions given. 
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4«2o4 Radiol sotope- heated thin shell system 

Without the solar heat flux# the above thin shell will 
gradually reach a very low uniform temperature. The solar 
Input not only heats up the shell but also create an 
asymmetry of temperature distribution so that a net 
reradlative force will act on the shell to generate a 
torque. If the shell is Internally heated to an elevated 
temperature# the superposition of the skew sinusoid (due to 
the solar input) on the elevated uniform temperature will 
increase the net reradiative force because of the higher 
energy of the reradlative photons. 

Practically# all alpha energy can be absorbed by a very 
thin layer of material# therefore alpha-emitting 
radioisotopes are good internal heat sources, Po-210# 
Pu-238# Cm-242 and Cm-244 are some of the possible 

candidates because they give the low radiation hazard. To 
analyze this design proposal numerically# a radioisotopic 
heat flux which corresponds to the amount of radioisotope 
deposited per unit shell area is added to Equation (4,2), 
The results of the optimum spin rates and the maximum 
erecting forces (normal component of the net force) as 
functions of radioisotopic heat flux are shown in Figure 
4,8, The addition of uniform heat flux favorably increases 
the optimum spin rate but unfortunately# the associated 
increase in the steady state erecting force is only about 
twenty percent as shown on Figure 4,8, (Note that the 
amount of isotope required varies with isotopes. For 
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Po-210, the typical isotope requirement using Table 5.1 for 
conversion is about a few mg/cm^of shell surface area.) 

4.5 Thin Shel 1 Approx Imat ion Method 

4.5 .1 Description of Thin Shell model 

in the Strip Approximation method/ both the radial and 
circumferential heat conductions are neglected. Since the 
circumferential heat conduction will tend to smooth out the 
asymmetric angular temperature distribution/ the Strip 
Approximation may overestimate the maximum and minimum 
temperature and hence the performance/ O.e./ the torque 
generated) « To account for circumferential heat conduction/ 
the governing differential equation becomes 

f C |X = K gr, + Jp [I, F( -- e «« 

where P/ C/ k, € , c~, lo /^t*/ and R are the density/ heat 
capacity/ conductivity/ emissivity/ Boltzmann constant/ 
solar constant/ shell thickness and radius 
respectively. F(B/t) Is the hal f *rectl f I ed sinusoid 
which characterizes the solar intensity around a 
cylinder; is the complement angle between the spin 
axis and sun~llne. 

This semi-linear partial differential equation is solved 
numerically and analytically in the following sections. 
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4«3a2 Semf 'irapllcl t numerical solution 

The Implicit difference formulation of the above 
partial differential equation requires linearization of the 
radiation term I . e . Th is non-linear 
term is not expected to Introduce serious numerical 
instability problems^ but the resultant coefficient matrix 
is not exactly tri diagonal because there are two non-zero 
terms on the extreme corners of the matrix due to the 
closed-loop nature of the circumferential heat conduct I ono 
As a result/ the matrix factorization technique cannot be 
applied and the Gaussian reduction method Is used instead/ 
which Is more time-consuming. The numerical results for the 
thin metallic shell using Thin Shell Approximation agreed to 
about one percent of the results of the Strip Approximation 
(see Figures 4.2/ 4.3/ and 4.4). 

Two of the by-products of this computation are the 
parallel and the normal force components (relative to the 
sun axis line) as functions of time. They both exhibit 
damped oscillation as shown in Figures 4.6 and 4.7 (for a 
shell length of a centimeter). From the spacecraft dynamics 
viewpoint/ the absolute magnitude of the normal force 
component Is Just as important as its magnitude relative to 
the parallel force component/ because the parallel and the 
normal forces are directly related to the precession torque 
and the erecting torque respectively. The following section 
discusses a proposed scheme to increase the normal force. 
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4»3o5 Subi iming material coated shei 1 system 

The thermal force on the shell originates from the 
recoil momentym of photons which have a spectrum of energy 
characterized by the temperature of the radiative surface. 
At steady state^ any unit surface on the shell will radiate 
as much energy as it receives during one spin cycle. The 
scalar momentum Is directly proportional to the energy 
radiated but inversely proportional the particle expulsion 
speed. Consequently^ if sublimation of a solid material is 
employed as the mode of removing energy^ the expulsion speed 
('^3*10 cm/sec) Is about a million times lower than the 
speed of radiated thermal energy (3*10'°cra/sec); the momentum 
resulting from sublimation Is therefore about a million 
times the thermal radiative momentum, Thus^ there is quite 
an incentive to coat the thin shell with subl iming material. 
The heat flux carried by sublimation can be derived (see 
Appendix L) as 

where is heat of sublimation 
wt i s molecul ar wei ght 
A/6 are the sublimation constants 
T Is the sublimation temperature 

which is added into the bracket of Equation (4,4) to 

account for sublimation. 

The resulting temperature distributions are shown in 
Figures 4.9^ 4.10 and 4,11, The important features are (1) 
the damped oscillations of the instantaneous maximum and 
minimum temperatures are essentially suppressed by the 
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sublimation heat loss for wide range of spin rate, however/ 
the instantaneous phase lag still retains a small damped 
oscillattcm; (2) the steady state temperature distribution 
Is lower than the results from a bare shell; (3) the phase 
Ifg Is not as large as in the case of a bare shell. For 
fast spin rate (oj>> cOopt the instantaneous normal and 
parallel force ccMuponents will tend to oscillate for a very 
long time because the surface tenperature does not have 
sufficient time during each cycle to equilibrate. There are 
two Interesting features worth mentioning: (1) for fast 
spin/ the restoring torque can actually switch polarities as 
shown In Figure 4.12; (2) as shown In Figures 4.13 (a and 
b}/ the amplitudes of the force components for low initial 
temperature will have growing oscillation for many cycles 
before the amplitudes of oscillation start to decrease. The 
Initial maximum temperature oscillates with damping/ but its 
mean value continues to rise toward the equilibrium maximum 
temperature. Since the net reradiative and subliming forces 
are proportional to the exponential power of temperature/ 
the amplitude of the oscillation attains a peak. 

4.3.4 Approximate analytical solution 

At steady state/ the angular temperature distribution 
for a bare thin shell as shoWn In Figure 4.14 has a 
remarkable resemblance to a sinusoid. This temperature 
profile travels around the shell with an angular speed of 
the spin. The Fourier expansion of this wave form for 
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constant spin (i) is 


1'Ct, 0) = 290'9+2/. 

-i2.i5 coszcu^ue:> -f i.^o-^iN^cujt+e) C4.&; 

+ 0'43'C0S36o;t+W 4 0-25 S/l\/3(Mlt-t-&)+-" 


which Indicates that the steady state temperature 
distribution can be approximated by 


TCt.6)= ISig £ ^ +A^, cosCU}t+9) + M2.S/NC0t+e>] 


= Tav3 [t + M 


C4.7; 


Here# is the average temperature# M is related to the 
maximum temperature difference and '^is the phase lag of the 
maximum temperature from the direction of maximum solar 
Intensity* The half rectified sinusoid F(€#t) and in 
Equation (4*4) are approximated by 


(-(tjd) 


2 ^ 

- / s/NC(Xt+e)-'^^ 

I \ 


C0S2iyi Ci^P+9^ 
4-n^ —I 


(T4.8J) 





Substituting the above expressions into Equation (4*4) and 
equating the coefficients of the constant# sine# and cosine 
terms# the solution of the resultant simultaneous equation 

T^% = l-'D/c'nof'^ C4.10) 

M = 0-SD /[ + C4 II) 

o|j = C4I^) 
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and the constants in Equations (4ol0}, (4/11) and (4.12) ai*e 
defined in Equations (4.17) through (4.22), 

For the radioisotope-deposited thin shell/ a constant 
heat flux term G is added to Equation (4,4) as proposed in 
Section 4.2, 4o The expressions for M and'f are the same as 
above/ however the expression for Is 


In the case of the subl iming material coated 
sub 11 mat I on heat flux term E 
to Equation (4,4), The above paratneters become 

cl^va +D/tr+lZy £■■ 10^ = 0 


shel 1 / a 
i s added 


C4, 14) 


M 


r= D/fe ^ 8 +4C + E ■ ( 


a 




2..5T 




nj' 


2 


C4/5-) 


^ f-Qfn 


-I 


f\ "TayQ ^ 


■eXj +‘^ c-e++f 


C4/6; 


The constant A/B/C/D/E/G are the coefficients of the 
partial differential equationo They are 


f\=yc C4.li) 

B = k/R^ C4-m 

C=-6crAr ^4,/9; 

D = JoCOs4^/^r C4.2Q) 

E = Hs my C 1144 ,14 2 1) 


97 



where % Is the heat of subll (nation 

RIhs is the radlolsotc^lcs heat flux 
and the rest are defined In Eq,{4,3) 

These approximate anal ytl cal expressions of the form of 
Equation (4»7) are capable of representing the true 
(numerical) tenperature distributions within about 20 % 
accuracy« In general^ Equation (4« 7) should be used with 
caution because it could introduce large errors in computing 
Other quantities such as the reradiative force components 
(note the reradiatlye force Is proportional to the 4th power 
of temperature). 


When the approximate analytic expression for 
temperature^ Equation is substituted In the normal 

force expression^ Equation (4.3)# the resultant expression 
Is a sixth order polynomial In terms of the spin rate a) . 
Applying the Oescartes Rule of Sign# it cari be shown that 
there Is one and only one real positive root In the 
solution. Analytically this root Is found by transforraing 
the sixth order pel ynomial to a cubic equation before using 


the Cardan * s f ormul a to ge t an expl i c i t express i on (HI) 

.^J/ 4 - r i -^1^1 /—I * 


The result Is 

2 


. . ^ _L_ f 2 c Li 

CUopt - i27f^ I 

+ SI f ^ q (f-f. 


cohere# 


% = U co^<l>/C2^r~R^3 > 
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For the bare and the. heated shells the nuraericat. results 
from this expression .asreed almost, exactly with the only 
positive root of the original ^Ixth order polynoraral. For a 
bare thin shelly the optimum spin Is nearly Inversely 
proportional to the shell thicicness and Is only a few 
hundredths of a radian per second as shown in Figure 4.15. 
For the radlolsotope‘:‘heated thin shelly the optimum spin Is 
greater by an order .of magnitude as shown In Figure 4,6. 
Agaln<> the approxtraate solution checked out well (X2S%) with 
the exact (numerical trlal-and-error) solution. 

For the case of the subliming material coated shelly 
the approximate optlmimi spin Is not obtainable because the 
subliming force is a complicated function of temperature. 
As a result, the integral of Iquatlon (4.3) is very 
difficult to obtain analytically. However/ based on the 
reradtative force results, the optimum spin is at least an 
order of raagnitude less than that of a bare thin shell. 
From a practical point of view, the coated shell functions 
best If the spin is about a few rpm because at this spin 
rate, the average parallel force Is quite small and yet the 
normal force for erecting the spacecraft is still 
appreciable as indicated in Figures 4.13 (a and b) . 

4.4 Thick Shell ADProxlmation Method 

For a thin shell/ the radial temperature variation is 
negltglble. However for a shell coated with subliming 
material/ the above assertion is not true anymore/ because 
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the thermal conductivity of the subliming material Is about 
a thousand times smaller than the metal and Its thickness 
has to be at least ten times as thick as the metal In order 
to be practicalo Fortunately^ due to the low thermal 
conductivity of the subl imlng material^ the circumferential 
heat conduction effect can still be Ignof^ as In the Strip 
Approximation method (see Section 4.2), The resultant 


governing differential equations and boundary conditions are 

£JL ^ is, £E, 

d ^ ar k at 


= - (T& 77^- Hs 
R+<ir 77/4- 


Bji;' 




j 77 Cn 0 ) = com 77^-a/t 


C4-.24) 


I ^ ^ 2 ^ 3 ^ d o. 6 d 

Since the shell Is divided Into U strips^ for a 
complete angular temperature distribution at any time t^ the 
above differential equation has to be solved <J tlmes.« 
Obviously^ this is a very 1 1 me- consigning processo As a 
matter of fact/ the number of operation steps will exceed 
the solution of the original two-dlmensIon parabolic partial 
differential equation (Equation 4«1} by ADI (Alternate 
Direction Implicit) scheme (see Section 5,4o3} for each time 
step because the number of radial mesh points Is generally 
less than the number of circumferential mesh points. The 
Important difference Is that the proposed method is very 
stable numerically. On the contrary/ the ADI method will 
likely introduce a severe numerical stability limitation due 
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to the non~reGtangu1ar geometry of this problem (j^) 6 As a 
result^ the present nt4mer I cal method will probably require 
less computer time for a complete transient solution in 
comparison with the ADi method because the present method 
can utilize a large time stepo 

There are many interesting features contained in the 
numerical solution as displayed in Figures 4.16 and 4.17o 

(1> The teroperature gradient across the shell is 
usually very significant (about 12 *‘C/cm) . 

(2) For each striP/ this gradient will switch signs for 
each spin cycle . 

(?) The phase lag of the maximum temperature on the 
subliming surface is less than the phase lag on the metallic 
surface by above 50|o 

(4) The temperature distribution at any instant 

resembles a sinusoid as shown in Figure 4.16, 

(5) The maximum and the minimum temperatures as 

functions of time shown In Figure 4.17 are similar in shape 
to the curves from the Thin Shell Approximation, (see 
Section 4.3,3} 

(6) Finally^ the computer execution time is about five 
times greater than that for the similar calculation by the 
Thin Sheil Approximation* 

4 .5 Comparison of Results 

The original two di mens ion time dependent heat flow 
equation has been solved with varying degrees of 

IjOI 


sophlstlcatlono it Is obvious that as each assumption Is 
removed/ the fine structures of the ’exact' solution will 
begin to unfold^ and the computer execution time will also 
Increaseo in general^ the Strip Approximation Is very 
accurate for thin metallic she1l« If the metallic shell Is 
not very thin# then the smoothenlng effect due to the 
circumferential heat conduction should be taken into account 
by the Thin Shell Approximation, When the shell exhibits 
poor thermal conduction such as with subliming material# the 
Thick Shell Approximation should be used In order to yield 
accurate surface temperature# because the sublimation 
process is very strongly dependent upon the surface 
temperature. Only under extreme conditions# such as a 
thick metallic shell coated with organic subliming material# 
must the original tvyo^dlmensl on differential equation be 
solved rigorously. Fortunately thl$ situation does not 
exist In this study. 
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FIG 4.2 TRANSIENT MAXIMUM 
TEMPERATURE OF BARE SHELL 
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FIG 4.3 TRANSIENT MINIMUM 
TEMPERATURE OF BARE SHELL 
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FIG 4.4 TRANSIENT MAXIMUM TEMPERATURE 
DIFFERENCE OF BARE SHELL 



106 


100 150 200 250 300 350 

TIME t, sec 



FIG 4.5 TRANSIENT PHASE 
LAG OF BARE SHELL 
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FIG 4.6 TRANSIENT NORMAL FORCE 
COMPONENT FOR BARE SHELL 
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FIG 4.7 TRANSIENT PARALLEL FORCE 
COMPONENT FOR BARE SHELL 
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FIG 4.8 

OPTIMUM SPIN AND MAXIMUM NORMAL 
TORQUE FOR A HEATED SHELL 
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FIG 4.9 TRANSIENT MAXIMUM 
TEMPERATURE OF COATED SHELL 
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FIG 4.10 TRANSIENT MINIMUM 

TEMPERATURE OF COATED SHELL 
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FIG 4.11 TRANSIENT PHASE 
LAG OF COATED SHELL 
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FIG 4.12 TRANSIENT PARALLEL FORCE 
COMPONENT FOR COATED SHELL 
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FIG 4.13a NORMAL FORCE COMPONENT 

FOR FAST SPIN ( w»o),,t ) 
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FIG 4.14 STEADY STATE TEMPERATURE 
DISTRIBUTION AROUND A BARE SHELL 
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FIG 4.15 OPTIMUM SPIN RATE DEPENDENCE 
ON BARE SHELL THICKNESS 
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FIG 4.16 OUASI-STEADY TEMPERATURE 
DISTRIBUTION AROUND COATED SHELL 
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te;;mpe:rature. 


FIG 4.17 TRANSIENT MAX AND MIN 
TEMPERATURE FOR THICK-COATED SHELL 
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CHAPTER 5 


performance analysis of microthruster 

5^1 Microthruster Design Considerations 

In conventional cheraical or nuclear rockets^ the 
thermal power from chemical reactions or nuclear fission 
reactions is Imparted to the propellant in the form of 
random gaseous kinetic energy which is transformed to 
directional kinetic energy by means of convergent-di vergent 
nozzle to produce thrusto Theoretically^ if all the gaseous 
molecules are expelled exactly In parallel with the thrust 
axis^ the resultant thrust will be the maximum attainable 
for a given chamber stagnation condition. 

In the subliming microthrus ter^ this maximiun can be 
achieved In principle If the subliming solid is placed at 
the focal point of a perfectly ref lecti ve parabol Ic nozzleo 
In reality^ this Is not a practical system because the 
subliming surface has to be very small at the focus and also 
the practical parabolic surface is far from being 

reflective. Nevertheless^ a simple system^ such as a 
cylindrical cup filled with subliming material as shown In 
Figure 5, lA^ can realize about half of the theoretical 
maximum thrust. The performance of this simple design can 
be easily improved by connecting a divergent conical nozzle 
to the cyl indr leal cup (Figure 5»1B), The plausible reason 
is that gas molecules sublime^ on the average^ with a speed 
of^-^:^] (^) which is higher than the sonic (pressure, wave) 
speed of(^^)(for complex molecules V«1,0)^ so that the 
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nozzle has to be divergent to increase the resultant thrusts 
Physically,; the molecular interaGtions with nozzle wall are 
conveniently divided into specular and diffuse reflect ions« 
In a divergent nozzle,; the wall is slanted outwards From a 
simple geometric consideration,; it can be shown that after 
either type of interaction with a divergent wall# the 
average final normal velocity component (parallel to the 
conic axis) is always greater than the Initial » 
Consequently# the longer the nozzle# the higher the thrust 
will beo However# the net mass flow# due to the nozzle 
obstruction# Is decreased to an asymptotic value as shown In 
Figure 5o2o To Increase the thrust further# the alternative 
is to Increase the mean speed of the molecules by heating 
the nozzle skirt so the molecules adsorbed by the wall will 
be evaporated with a higher average energy than their 
initial average energy o 

5.2 Radioisotope ApdI icatlons 

5(,26l Radioisotope properties 

There are at least two areas where radioisotopes can 
contribute. These are as the heat source for the 
sublimation process and as propellant for thrust generation 
by the deeay particle recoil mpmentumo Unlike Joule heating 
(electric resistance heating}# the radioisotopic heating 
cannot be turned on and off so that It imposes a design 
problem^ On the other hand# it is very reliable in 
generating a predetermined thermal energyc Considering 
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decay recoil momentum^ the heavy particle eniltter Is very 
attractive for generating small thrust as will be shown In 
Section 5o2o2« In the present study^ only alpha-emitters 
appear to be useful because of both their ease of converting 
the decay energy to heat and their avallabi lltyo The 
selection of a particular alpha-eral tter depends on many 
factors such as half-life^ radiation hazard^ cost^ and 
mater I a1 compat I b! 1 1 1 yo 

In Table a few radioisotopes of interest are 

listed along with their physical properties« These Isotopes 
are formed by capture of one or more neutrons o The nuclear 
rad iatlons given off by them emanate directly from their 
natural decay scheme or from Interaction of primary 
radiation with other materlalSo Usually^ only gamma rays 
and neutrons are Important In the required shielding 
considerations since they are more penetrating than both 
beta and alpha particleso Although the isotopes of 
lns>ortance In this study decay by emission^ gamma rays and 
neutrons are present due to impurities^ ( ^n) threshold 
reactions with tow atomic weight elements (e6go^ 0-17 and 
0-18)^ spontaneous fission (except Po-210)^ and induced 
fission Ceogo^ Pu-238 and Cm-244)o The effects of neutron 
multiplication (jL^) In Pu-258 and Cm-244 are negl igl bl e for 
the present applications^ however^, the effects of impurities 
are sometimes unavoldable» Incidentally^ Th-228 is an 
excellant propulsion Isotope; the half-life is 1«9 year and 
Its decay causes emission of five energetic alpha particleSo 
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Unfortunately^ It Is not on the AEC isotope production list. 

5.2.2 Radiolsotope°^doping for fast despin 

5ubliroat!on processes^ as with other phase transition^ 
re qui re thermal energy. For materials of interest^ the 
sublimation heat Is of the order of 15 Kcal/mole 
Unfortunately^ subliming material are usually very poor 
thermal conductors^ so that to supply a large amount of heat 
by way of conduction through these materials is very 
undesirable because of the associated large temperature 
gradients and the non°° uni form subl iming surface temperature. 
Radiol sotope^doping seems to offer an excellent solution. 
Minute amount of alpha°°emitting radioisotopes could be mixed 
in with the subliming material to form an internal heat 
source so that the sublimation heat is readily available 
and the surface temperature can be made fairly uniform 
radially. The resultant thrust would be about a few dynes 
per square^cent Imeter (based on the vapor pressure of Figure 
6. 5 at “275*K)^ which is sufficiently large to rapidly 
despin an initially fast spinning spacecraft. 

5.2.3 Radioisotopic decay for long term spin-control 

If the microthruster cup bottom is coated with a thin 
layer of an alpha-emitting radioisotope, then a small thrust 
will act on the cup due to the alpha recoil force as soon as 
the subliming material is depleted. The purpose of this 
film of radioisotope is to generate a long-term despin 


124 



thrust for the correction of any unexpected spin-drift* 
Since the alpha particles have such a low penetration 
capability^ half of them are absorbed by the bottcxn of the 
cup and the other half are assumed to escape through the 
flow channel* The thrust non-obstruction factor for these 
high energy particles Is very much different from that of 
the thermal energy molecules (see Section 3*2*2} because 
IntefaGtl ons with channel surface by these high energy 
particles will heat up the surface and the subsequent 
reemi sston will account for only a very small fraction of 
their Initial momentum (less than a thousandth of the 
original}* These high-energy charged particle Interactions 
with the channel surface will contribute negligible recoil 
force* Consequently^ as a good approximation^ the thrust 
for high energy particles can be assumed to come from the 
direct streaming alone (those experiencing no channel 
surface obstruction)* The thrust non-obstruction factor for 
high-energy charged particles (see Figure 5*3) Is a function 
of both nozzle length and nozzle cone angle and is always 
less than that for the thermal energy subliming vapor* 

The absolute thrust due to alpha particle emission is 
about three to four orders of magnitude smaller than the 
sublimation thrust* This is obtained in the following 
manner* The relativistic expression for the recoil momentun 
of a alpha decay Is 
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where KE is the alpha kinetic energy 
C is the speed of light 

Is the rest mass of alpha particle 

The total number of alpha emissions per second is equivalent 

to the activity of the radioisotope and it is given by 

dN/dt=-ANe-« 

H/fiere /V IS ^/cm^ 

Taking into account the factors of 1/2 for the 
absorption loss and 1/2 for surface normal component 
correction (not all particles leave the surface Vertically)^ 
a useful expression for a thin“layer nuclear (neglecting 
self absorption) disintegration recoil force Is 

Ft^) - o . 8 is3) 

where WPA is watts of radioisotope deposited 
on a em^ of surface V 

For WPA“0ol W/em^ ^ the recoil force is about 4o7*10 
dyne /cm*^ » With inclusion of the nozzle skirt this force 
will be reduced by a factor equal to the thrust 
o©n~obstruct ion factor^ 

However^ as the thickness of the Isotope layer (or WPA) 
increases further., the above expression will overestimate 
the thrust because of the self^shieldingo in principle^ the 
thrust will reach a plateau for isotope thickness beyond Its 
alpha rangOo if the isotope thickness is lese than the 
range the thrust Is given by (iJ|) 

1/2, 

T=(g|-) C-f'-W 
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For Po“210^ the graph of thrust versus the dlmensionVess 
thickness x/rg. Is shown In Figure 5o4o It Is evident that 
the self absorption limits the maximum thrust intensity to 

—3 

about 5*10 dyne/ cm*’ o In other words^ the thickness of 
deposited layer should not exceed the alpha range of about 
lo7*lQ^ cm 6 

For available radioisotopes such as Po-llO^ Cm-242^ 
Cm-244 and Pu-238^ the life span of their decay thrust, 
which is proportional to the half’^lives, vary from about a 
year to several hundreds years, so that this type of nuclear 
disintegration thrust is suitable for long-term spin 

control 0 

5o 264 Radioisotopic heating of nozzle wall 

For a given mass flow rate, in order to increase the 
thrust and hence the specific impulse, more energy has to be 
Imparted to the vapor « In an advanced nuclear rocket, 

seeded propellant can absorb thermal radiation from the 
divergent nozzle to increase the exhaust velocity, but here 
the sublimed vapor is so rarefied that the molecules do not 
interaet with themselves or with thermal radiation^ 
However, if the nozzle wall Is heated by radioisotope decay 
to a temperature much higher than the sublimation 
temperature, molecules adsorbed by the surface will be 
reemi tted with a higher average speedo Only the adsorbed 
molecules will benefit from the heated wall, but for a 
nozzle of moderate length the majority of the molecules will 
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Indeed Interact v\/ltfi nozzle surface^ The results of the 
heated~wall effect is shcvn in Figure 5=5^ for the case 
where the wall temperature is maintained at twice the 
subliming temperature« For the unheated-wal I case^ the wall 
temperature was assiMed to be the same the subl I tiding 

temperature^ but in reality the wall temperature would be 
lower so that the relative effect of having a heated wall Is 
actually higher than that shown In Figure SoS, The wall 
temperature Is dependent upon amount of radioisotope 
deposited^ surface emissivity^ radiation view factor and the 
heat loss to the reemitted moleculeso 

5 . 5, Xv 1 i ndr l ea I -cup M ic rotkruster 

5»3ol Description of governing equation 


The schematic diagram of this device Is shown is Figure 
SolAo The sublimation heat is supplied from the cup bottom 
where the heat is generated by a few layers of radipi sotopeo 
The outside cup surface Is assumed to be therraally insulated 
CioCo^ negligible surface emissivitylo Since sublimation 
Involves both the heat transfer and the mass transfer^ the 
temperature distribution is always In transient stateo The 
subliming surface temperature^ and hence the resultant 
thrust^ as a funetlon of time must be obtained from the 
solution of time a dependent one^dl mens ion heat flow 
equationo The governing equations are the followings 


3 ^* 


;>cil 

^ at 

dX,/dt= - GCi)/P 
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The ass0c Sated boundary conditions are 

-|f£E| = ^cti-€,TT^coA) 

9X, IX’*o 

CS".' 

XoCO )= constant 

where GCT) Is the sublimation rate 

q(t) Is the radioisotope heat flux 
fCx) Is the Initial temperature profile 
Tfi Cfi Ha are the density^ heat Gapaclty^ 
thermal Gonduct I vl ty^ and heat of 
sublimation respectively 

Two numeriGal approaches^ the explicit and the Implicit 
difference methods^ have been trled<, 

5o3o2 Explicit solution of temperature distribution 

The difference formulation Is such that the future 
temperature at any point Is simply dependent upon the past 
temperature only; however^ the numerical stability criteria 
requires the time step and the mesh size to be chosen In 
such a way that AtMx^ Is always less than or equal to 
(Oo5}(^c5/k (about 200 In cgs unlt)o This limitation 
complicates the problem somewhat^ because physically the 
outermost mesh will always decrease In size due to 
continuous sublimation process o This difficulty Is resolved 
by starting with a larger than minimum outermost mesh; then 
as soon as the outermost mesh reaches the minumum allowed 
size^ It is combined with Its neighboring mesh to form a new 
outermost mesho 
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Th© result ©f computatfon for one design Is shown in 
Figures 5 <,6^ 5o7^ and 5o8o The initial temperature was 
assumed to be uniform (not true In reality)^ The surface 
tempera ture dipped to a minimum at the beginning because the 
heat losses through sublimation and radiation are larger 
than the heat supply through conduction (almost zero) and 
then started to rise to a quasi °°equl 1 1 br I um value because 
the resultant large temperature gradient due to Initial dip 
Increased the conduction heat transfero The surface 
temperature crept up slowly to a maximum as a result of the 
reducing thermal resistance to the flow of heat from the 
radioisotope layer to the surfaceo An soon as all subliming 
material was evaporate# the teraperature of the bottom of the 
cup rose quickly as shown In Figures 5a9# SolO# and Soil# 
because the only heat removal mode was then thermal 
radlatlone, However# an elevated equilibrium temperature was 
rapi dly reached, 

5<,3o3 Implicit solution of temperature distribution 

To avoid the restriction on small time steps and the 
coarse mesh size# the implicit difference formulation Is 
employed. However# the temperature at a point depends upon 
both the past temperatures around this point and the current 
temperatures around this point which are unknown. 
Consequently# the spatial temperature distribution at any 
time Is the solution of simultaneous linear algebraic 
equations. Fortunately# the coefficient matrix of these 
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simultaneous equations is tridiagonal so that the very 
efficient matrix factor I zatlon teGhnique can be emp1oyed« 
The solution is very stable In general o Externally 

Introduced disturbances such as the rearrangement of 
boundary and mesh net are usually damped out in a fraction 
of a minute (real tilie) o The results of a sample 

calculation based on the same conditions as In the case of 
the explicit solution are shown in Figures 5,9^ 5ol0^ and 

Sill for the time period during which the subliming material 
is about to disappear completelyo 

5«3o4 Conparlson of results 

The results of the above two methods agree very well 
throughout the tested time span except at a few transition 
instants when the mesh net was rearrangedo As long as the 
surface was coated with subliming material# the surface 
temperature remained relatively constant after the initial 
fast transient because of the large damping effect of the 
sublimation process (see Figures 5o7 and 5ol0)o Since the 
thrust increases exponential ly with the surface temperature# 
a small variation in temperature can be important. 
Consequently# accurate prediction of the temperature is a 
prerequisite for accurate evaluetlon of the thrust. Between 
these two numerical methods the implicit scheme is preferred 
on this count. An interesting fact found in testing the 
above two programs was that the explicit computer coding 
required many checking calculations to prevent instability 
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and as a result actually requi red more computer time than 
the In^llclt program under more or less the same condltlonSo 

5»3»5 Thrust non-obstruct Ion factor 

By analogy with the mass non-obstructlon factor^ the 
thrust non-obstruction factor Is the ratio of real thrust 
from a channel to the thrust from a plain subliming surfaceo 
Unlike the mass non-obstruction factor# the thrust 
non-obstruction factor can be# and Is desirably# greater 
than unityc These tw non-obstruction factors are different 
because the thrust factor Involves a (v*cos3) weighting 
factor In the integration of the Maxwe.11ian di strlbutIon« 
Physically# molecules specularly reflected from a vertical 
wall do not affect the normal force componento However# 
those adsorbed and subsequently reemltted diffusely will 
lose some Initial normal momentum on the average# so that 
the thrust non-obstruction factor for a cylindrical channel 
will be less than unity In general (except when the 
evaporation coefficient Is very small )<, Also# since the 
average wal l-!ncident angle (with respect to channel axis) 
Is larger than the average angle for molecule escape from 
the channel wall# the thrust non-obstruction factor Is 
expected to be larger than the corresponding mass 

non-obstruction factor » These effects are shown In Figures 
5ol2 and 5ol3 where non-obstruction factors are plotted as 
functions of the L/R rati© and the evaporation coefficient# 
(forO(®loO# the sublimation surface specular reflection Is 
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2ero>6 It Is Interesting to note that for small evaporation 
coefficients (large subliming surfaee reflection) the thrust 
non-obstruction factor can actually be greater than unity 
due to the large subliming surface reflection, 

GhoBPed-CQne MI croth rus ter 
5,4,1 Introduction 

To significantly increase the thrust non-obstruction 
factor^ the flow channel (or nozzle) has to be modified. 
Intuitively^ a divergent nozzle should Increase both the 
mass flow and thrust because the slanted channel wall helps 
to reflect molecules outward. Considering fluid dynamics^ 
there Is also a plausible reason. Vapor molecules leave the 
subliming surface with an average speed of (1TkT/2m) which 
Is greater than the sonic or pressure wave speed Qkl/m) of 
any non-monoatomic gas (noting the specific heat ratio 
»l+2/f, where f Is the number of degrees of freedom) (51) . 
The flow speed In this supersonic region can only be 
Increased If it Is expanded through a divergent nozzle. Of 
course the equations developed In fluid dynamics cannot be 
applied directly to the present pressure range^ but 
apparently.^ useful qualitative conclusions can still be 
drawn from them. If the above deduction Is correct then.^ a 
thrust engine design Involving the heating a subliming 
material In a coated chamber to produce vapor to be passed 
through a nozzle (JJ.) would not result In an efficient 
desi gn. 
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The simplest design of a divergent nozzle Is probably a 
chopped coneo However, It will be shown that proper design 
can Improve Its perforreanee significantly such that the 
thrust does actually approach the theoretical maxlroumo Two 
of the most Important geometric design parameters are the 
optimum cone angle and the optimum nozzle skirt lengtho In 
theory, the longer the skirt length the higher the exhaust 
speed (or the better the propellant expansion), however, 
"diminishing returns'* will limit the lengtho In Figure 5,2, 
It Is obvious that contributions to the thrust are 
decreasing for L/R greater than about six for a 45-cone 
nozzle« A schematic diagram of a cylindrical-cup 

microthruster is shown In Figure SolB,, 

5o4o2 Optimum cone angle calculation 

The optimum cone angle Is very difficult to obtain 
analytically from Clauslng's integral equation (see Section 
3o2) The main difficulty lies In the formulation of the 
probability kernels, especially for slightly complicated 
assumptions on gas-surface Interact lonso However, 

numerically the non-obstruction factors for a chopped cone 
geometry pose no additional problems other than those 
described In Chapter 3, In the cylindrical tube case, the 
point sources from the Inlet (apparent throat) were divided 
Into differential cones, which Intersect the channel wall or 
trace space curves on the wallo The number of molecules 
contained in each cone are proportional to the spherical 
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cosine distribution so that the surface adsorption density 
In each channel ring is the properly weighted sumation of 
ail intercepted space curves from all cones of the point 
sources o 

Unfortunately^ truncation errors introduce an 
uncertainty by at least a fraction of a percento The major 
portion of the uncertainty probably comes from the numerical 
treatment of the Interceptions of differential cones by each 
channel rlnga The computer program (CONNOF) for solving the 
Ghopped cone geometry treats the differential cones in a 
slightly different fashions from the cylinder case (TUBNOF) 
so that a direct ccxnparison can be made with the 
cyl indri cal “tube solution (noting a 90 degree chopped cone 
is equivalent to a cylinder)o 

In this computer program (CONNOF)^ the circumference 
of each individual differential cone is equally divided 
into about one hundred equal groups such that each group is 
considered as a beam bundle characterized by a set of polar 
and azimuthal angleso The trajectory of this beam bundle Is 
traced until it Is adsorbed by the channel wal 1 « Adopting 
the same procedure as In the cylindrical tube case (see 
Section 3o2o2)^ the mass non"Obstruct Ion factor is conputedo 
Since the polar angle of each escaping beam is known^ the 
normal momenti^ and hence the thrust non“Obstructlon factor 
can also be computed« For a given nozzle lengthy the 
non-obst ruction factors are calculated for different cone 
angl es o 
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The results are shown in Figure 5ol4o For maximum 
thrust there Is an optimum cone angle which Is strongly 
dependent on the gas-surface properties^ in particular^ the 
specular ref lecti v! ty« In Figure 5 <,14^ although the thrust 
reaches a peak^ the mass flow monoton i cal ly increases as the 
cone angle Is widenedo Also In comparison with the results 
from TUBNOF^ it Is surprising that the two results agreed 
very well (see Table 5<.2)o By manipulating the cone angle 
of a short nozzle the thrust level can be increased by about 
twenty-five percent as shown in Figure 5ol4o In Figure 5 ^ 2 , 
It is shown that by lengthening the nozzle^ the thrust can 
be increased by about twenty percent « Consequent! it 

pays to design a chopped cone microthruster properly because 
it Is capable of Improving the efficiency Gver the 

cyl Indrlcal-cup design by about fifty pereento With a 

radioisotope-heated nozzle skirt as described In Section 
the performance can be Increased further by about 

another twenty-five percent as shown In Figure 5o5o 

5«4o3 ADI solution of temperature distribution 

In Section 5<,3^ the temperature distribution In a 

cylindrical cup was obtained from the solution of the one 
dimension parabolic partial differential equation by 

assuming negligible radial heat transfer. The same 

assumption is not reasonable for a chopped cone geometry^ 
because the radioisotope heat flux will funnel out from the 
bottom and introduce radial heat flux component even with 
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perfect thermal insulation of the external cone surfaceo 
Consequently^ the temperature distribution in a chopped cone 
Biicrothruster must Include the effect of radial heat 
conductlono A straight Implicit difference scheme would 
require an Inversion of about 100*100 matrix for each time 
stepo Thi s method Is certainly not practical because of the 
enormous computer time needed as well as the large but 
unknown numerical error associated with this type of 
operation (Inversion of large matrix) o A relatively new 
technique called ADI (A1 ternate-Dl rect lon-lmpl i ci t) was 
developed for fast convergence (50) of parabolic partial 
differential equatlonso This technique makes full use of 
the fast matrix factorization operation for solving the 
resultant matrix equation and is unconditionally stable for 
any mesh size or time step provided rectilinear coordinates 
are used for the partial differential equatlono However^ 


the geometric symmetry of 

the 

chopped cone 

demands use 

of 

cylindrical coordianteso 

The 

ADI method 

was chosen 

for 

solution using cylindrical 

coordinates 

wi th 

the 


understanding that the method might not be uncondl ti onal 1 y 
Stable In this coordinate systemo The description and 
listing of the computer code are given In Appendices Jo3 and 
Po7» 

While testing the computer program^ the implicit 
solution of the center-line in the axial direction exhibited 
an I r regul ar i ty whose cause has not been determined^ but by 
a parabolic interpolation of the center-line temperature^ 
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this difficulty has been bypassed^ As suspected^ the 
numerical scheme Is not unconditionally stablso For a 
reasonable fine mesh net (0<,l*0ol crn^ the corresponding 
max Imum time step should be less than about 1.0 second. As 
a result^ the computer execution time on IBM 360/65 Is 
comparable to the real time (total subliming time} thus for 
a propellant Csubliming material) thickness of a few 
centimeter^ several hours of computer time will be required 
to complete the transient temperature solution. The typical 
temperature distribution shown In Figure 5.15 suggests that 
the one-dlmensIon solutions of Section C5.3) may be used as 
a first approximation because the temperature varies very 
little radially except In the ylclnlty of boundary. Only 
for the detailed design calculation will the computer 

program need to be executed In Its entirety. 

5.4.4 Performance parameters 

In general^ the thrust level for this type of 

microthruster Is less than a hundred dynes or a few hundred 
micro-pounds. For ease of comparison with other thruster 
systems^ the parameters characteristic of a thruster were 
calculated. 

Neglecting the pressure differential effect at the 
exit# the exit velocity Is just the ratio of thrust to mass 
flow rate. The exit velocity as a function of the nozzle 
length Is plotted In Figure 5.16. Since the exit velocity 
Is supersonic throughout the nozzle length# the throat 


138 



cross”sectf on Is not a real quantity^ but Its effective 
location can be found by extrapolating the exit velocity 
until the sonic velocity Is reached as shown by dashed line 
In Figure 5ol6<, The corresponding fictlous throat area At f 


the thrust coefficient and the specific impulse Igp # are 
(naphthalene) 

At = 'h' ( R - i /'t'Q n ^ O'S A 0 S ) 

C = x/CB'At) = 1.8 C5:9} 

Isp= Ve/9o=21, C5>10) 


Since the performance Is strongly dependent upon the 
propellant molecular weight and the temperature range over 
which the ralcrothruster can be properly operated^ the 
selection of a suitable subliming material Is Inportanto In 
this study only the biphenyl and naphthalene are considered 
becau se of their adequate vapor pressures In the vicinity of 
273°K and their relatively high melting points. Some of 
their physical properties are given in Table 5.3o 
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FIG 5.1 SCHEMATIC DIAGRAM OF 
SUBLIMING MICROTHRUSTERS 



A 




B. Chopped Cone 


NON-OBSTRUCTION FACTORS 
vs L/R FOR 45* CONE 

'Moiovj Noiionyisao-NON isnaHi 


>1 'yoiovj NOiionyisao-NON ssvn 


LENGTH TO RADIUS RATIO, 


NON-OBSTRUCTION FACTOR, 


FIG 5.3 

NON-OBSTRUCTION FACTOR OF 
60“ CONE FOR HIGH ENERGY a 
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THRUST, 10" dyne/o.,2 


FIG 5.4 THRUST INTENSITY FOR 
PURE Fb“° DEPOSIT 
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THRUST NON-OBSTRUCTION FACTOR 


FIG 5.5 

HEATED NOZZLE EFFECT ON 
THRUST NON-OBSTRUCTION FACTOR 
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FUEL THICKNESS, 


FIG 5.6 REGRESSION OF FUEL IN 
CYLINDRICAL MICROTHRUSTER ( BOL ) 



TIME, sec 
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SURFACE TEMPERATURE, 


FIG 5.7 SURFACE TEMPERATURE OF 
CYLINDRICAL MICROTHRUSTER (BOL) 
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MAX TEMPERATURE, 


FIG 5.8 MAXIMUM TEMPERATURE IN 
CYLINDRICAL MICROTHRUSTER ( BOL) 
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FIG 5.9 REGRESSION OF FUEL IN 
CYLINDRICAL MICROTHRUSTER (EOL) 
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FIG 5.10 SURFACE TEMPERATURE OF 
CYLINDRICAL MICROTHRUSTER ( EOL ) 
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FIG 5.11 MAXIMUM TEMPERATURE IN 
CYLINDRICAL MICROTHRUSTER ( EOL) 
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NON-OBSTRUCTION FACTOR 


FIG 5.12 

NON-OBSTRUCTION FACTORS vs 
L/R FOR CYLINDRICAL CONE 



L/R RATIO 
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FIG 5.14 OPTIMIZATION OF CONE ANGLES 
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CONE ANGLE , 



FIG 5.15 SURFACE 
TEMPERATURE BY ADI METHOD 
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FIG 5.16 

EXTRAPOLATION OF THROAT AREA 
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TABLE Sol 


CHARACTERISTICS OF RADIOISOTOPIC HEAT SOURCES 



POLONIUM 

PLUTONIUM 

CURIUM 

CURI UM 


PO-210 

PU-238 

CM-242 

CM-244 

HALF-LIFE (YR) 





ALPHA DECAY 

0o38 

86 

0o45 

18 

FISSION 

00 

4o9*10'® 

7 0 2*10^ 

lo4*10‘ 

DECAY (MEV) 





ALPHA 

5o3 

5o49 

6 oil 

5,80 

GAMMA 

0o8 

0o044 

0,04 

0,04 

FUEL FORM 

METAL 

OXIDE 

0X1 DE 

OXIDE 

DENSITY (G/CC) 

9o9 

10 

9 

9 

W/GM (PURE) 

141o 

0o56 

120 

2,65 

PURITY 

0o95 

0o80 

0,90 

0,95 

W/CC (FUEL) 

815 

3o9 

882 

20,4 

Cl/W 

32 

30 

28 

29 

SHIELDING 

MINOR 

MINOR 

MINOR 

MODERATE 

CAPACITY OCG/YR) 

1 

15 

0,15 

3 

SOURCE STRENGTH 





(N/SEC-W) 

3o2*10^ 

5 0 2*10^ 

4,0*10^ 

4, 6*10"^ 

FUTURE COST 





($/W) 

20 

540 

17 

64 

•REFERENCES (54) 

. (li) 
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TABLE 5o2 


eOMPARlSON OF RESULTS FROM TUBNOF AND CONNOF 


L/R 

r 

Rffl 



(REF) 

(TUBNOF) 

(CONNOF) 

0o5 

0.0 

0.80291 

0.8023 

loG 

0.0 

0.67293 

0.6728 

1.0 

0.1 

0.69910 

0.6988 

1,0 

0.2 

0.72630 

0.7261 

2.0 

0.0 

0.51486 

0.5150 

0 

0 

0.0 

0.41974 

0.4211 

4o0 

0.0 

0.35619 

0.3589 

5.0 

0.0 

0.31074 

0.3138 


♦NOTE THE CONE ANGLE FOR CONNOF CALCULATION 


IS FOR 880 99 DEGREES RATHER THAN 90 DEGREES 



TABLE 5,5 


PHYSICAL PROPERTIES OF NAPHTHALENE AND BIPHENYL 



NAPHTHALENE 

Bl PHENYL 

FORMULAR 

C10-H8 

C12-H10 

M,Po (t) 

80,3 

00 

o 

00 

SPEC, HT, 



(CAL/GM-t) 

0o315 

0,307 

THo COND, 



(W/CM-t) 

0,0058 

0,0037 

MOLE, WT, 

128,2 

154,2 

HT, OF SUB, 



(KCAL/MOLE) 

14,23 

17,37 

DENSITY 



(6M/CC) 

1,145 

00 

o 

CONST, A 

11,45 

10,38 

CONST, B 

3729 

3799 

SONIC SPEED 



AT O^C (CM/SEC) 

1,367*10'^ 

1,243*10"’' 


•NOTE THE ABOVE CONSTANTS WERE OBTAINED PARTLY FROM 
PHYSICAL-CHEMICAL CONSTANTS OF PURE ORGANIC COMPOUNDS^ 
VOL II, 
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CHAPTER 6 


APPLICATIONS TO SUNBLAZER-TYPE SQ^AR PROBE 

6.1 Spacecraft Dynamics 

6.1.1 Passive attitude control system 

A complete passive attitude control system for the 
present Sunblazer probe would consist of an erecting shelly 
a restoring vane^ an annular viscous damper and two pairs of 
ralcrothrusterso The reason to have two pairs Instead of one 
pair Is mainly to maintain the ax I symmetry of the moment of 
Inertia. 

However^ a requirement for both the erecting shell and 
the restoring vane may be questionable. In principle^ 
either the erecting shell or the restoring vane should be 
able to provide the necessary eredting and precession 
torques to stabilize the spacecraft. The net reradlative 
force acting on the erecting shell Is shown In Chapter 4 to 
have a normal component and a parallel component. The 
normal component Is capable of aligning the spacecraft with 
the sun line such that the solar panel at one end of the 
probe will face the sun. The parallel component will cause 
the spacecraft to precess about the sun°*pointing axis. Even 
If the initial precession angle (angle between the 
sun**poInting vector and the spin vector) Is greater than 90^ 
the normal force component will/ due to the erecting shell/ 
Increasingly reduce the cone angle to less than 90*^. Before 
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the cone angle 

Is reduced 

below 

90* ^ the 

precession is 

retrograde and 

thereafter^ 

the 

precess i on 

becomes 

direct 

(p recess in the 

same sense 

as 

the spin)o 

This 

di rect 


precession Is similar to a spinning top under a gravity 
torque (l£)o As long as the precession torque does not 
exceed a Hmlt (see Section 6olo3)^ the spacecraft will be 
very stableo In the meantime^ the erecting torque^ though 
dirnlnlshlng in magnitude^ will gradually reduce the cone 
angle toward zeroo 

On the other hand^ a restoring vane can be designed to 
accomplish the same tasko This is done by making the 
external surface of the vahe^ which will "see" the sun# a 
moderate absorber of solar radiation and a poor thermal 
radiation emittero The internal surface^ however^ is coated 
with a highly emissive material so that the absorbed energy 
is radiated away from the sun by the Internal surface of the 
vane as shown In Figure 6olo In this manner^ the solar 
pressure on the external surface of the vane^ due to the 
absorption and the reflection^ will always sustain a stable 
retrograde precession because of the negative precession 
torque (see Figure 6o2)o The thermal reradiation from the 
internal surface of the vane will have the same function as 
the erecting shellp The net reradiative force has an normal 
component to erect the spacecraft and a parallel component 
to cancel only a portion of the solar pressure force acting 
on the external surfaceo 
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As a result^ both the erecting shell and the restoring 
vane are able to erect the spacecrafto The main differences 
between them ares Cl) the erecting force from the erecting 
shell will always be larger than that from the restoring 
vane on per unit surface area basis; (2) The erecting shell 
will force the spacecraft to have direct precession rather 
than retrograde precession as in the case of using restoring 
vanoo 


6«lo2 Euier^s dynamical equations 

The extension of Newton’s law to rotational motion is 
that the time rate of change of angular momentum is equal to 
the applied moment^ i«eo 




Ci-O 


where XL is the angular velocity for axis system. 

In an axisymmetric spacecraft# Ij^iy^lT (subscript T for 
transverse axis) and the relation between the absolute 
angular velocity of the body to# and that of the body axls#^ 
is 

tOy J i ^ 

The corresponding component form of the above moment’s 
equation# or the Euler’s dynamical equations# are 

Mx — Ix 

= ly d lO^/dt — SLx. -i-SLzIx 2^ 

|V1g — If 
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The uj^s are not angular velocity components in the ordinary 
sense of being derivatives of some spatial coordinates which 
describe the position of the body at time to To describe 
the orientation of the spacecraft^ Euler’s angles are used 
as shown In Figure 6o2o in terms of Euler’s angte^ the 
Euler’s dynamical equations become (37) 

= Ixdu;x/c/t 

= IrB -Irf ^ll^SCOSe+hO^^SINB 

M2- = IjdC^SiNB)/dt -0/x(>a.+elTf CO50 

• • 

tuhtte =: <p-t- ^ C05 B 

In principle^ the spin moment ^ the precession 

moment and the erecting moment M^-as a functions of time 
can be founds so that the orientation of the spacecraft at 
any time t is the solution of the above nonlinear coupled 
differential equationso The difficulties are; (1) the 
large amount of computer time required to obtain the moment 
components as a function of times (2) too many initial 
conditions invoivedo Consequently^ the numerical 

integration solution will be very t Irae°°consuming and the 
result will be difficult to analyse and to draw useful 
conclusions from o The following sections deal with only the 
specialized solutions to Euler’s dynamical equatlono 

6ol«3 Effect of various moments on precession 

For steady precession as shown in Figure 6o2^ 6 and wu 
are constant^ so the the precession is given by 


C6.3J 
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4 = i^±rr-^^f- 

^ ZljCOSQ i^ZIjCOSB^ IrSlNe 


l/z 


C6.4j 


For a given 8 and the precession* torque plot is 
shown in Figure 6«3o There is a maximum aiiowable 
precession (restoring) torque mJ ^ above which steady 
precession is impossible because the spacecraft will flip 
upside down^ where 

^ i~Qn 9 C9>5) 

The positive precession torque will force the spin axis 
to process about the preferred axis 1 (i»eo^ the 
sun*polnting axis) in the same sense as the spin^ as long as 
the torque Is less than The net force acting on the 
spinning spacecraft due to temperature asymmetry on the 
erecting shell will always have a positive precession torque 
componento Thus a solar vane or solar sail might have to be 
added to insure a stable precession^ loCo^ to balance out 
ail or part of the positive precession torquco 

In addition to the positive precession torque 
component^ a highly desirable erecting torque is also 
presents This erecting torque will force the spin axis to 
spiral into the preferred axis or the axis of precession^ in 
other words^ the precession cone angle Is continuously 
dlmlnishinge However^ unless the restoring torque decreases 
as fast as sin8^" a possi bll i ty exists that the precession 
becomes unstablee 
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Also^ during the stabilization period^ the spin rate 
slows down<j> which has two serious inpl icat Ions, First of 
a11^ the precession torque will increase with decreasing 
spin because the parallel force component of the thermal 
radiation reaches a maximum at zero spin. Secondly^ as 
shown in Equation (6,5) the maximum allowable precession 
torque decreases In direct proportion to the square of the 
spin. As a result^ to maintain a stable precession^ the 
despin ralcrothruster has to be designed Garefully^ l,e,# It 
should not despln the spacecraft rapidly to zero, 

6,1,4 Despln rate 

The spin behavior of an axi symmetric spacecraft is 
essentially dependent upon the microthruster spin moment 
alone. Its governing differential equation Is simply 

Ms = L.cJc^x/rft 

For an estimate^ assume there are two pairs of (i,e,# a 
total of four) raicrothrusters^ each of which has a throat 
radius of a centimeter^ the microthruster has a momen t arm 
of 20 cm and Its surface Is maintained at about 273’ K, The 
despln moment Is about 

M«“4*20 cm*15 dyne=-1200 dyne-cm C(>-7) 

A typical small satellite has an axial (or largest) moment 
of Inertia of about 

lx *5*10^ gra^cra C^-8) 

Consequently^ the despln rate is of the order of 2,4*10^ 
rad/ sec which is large enough to despln the Sunblazer 
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satellite froia 200 rptn to few rpm In 20 hours without the 
help of a Yo-Yo deviceo However^ about a pound of subliming 
matertal (obtained by multiplying the total subliming 
surface area<^ surface regression rate^ despin time and 
material density) is required to accomplish the task/ so 
that integration with a Yo-Yo is needed to reduce the fuel 
requirement to a tenth of a poundo As to the residual spin 
torque from the alpha emission at the bottom of the cup 
after the material had sublimed away^ for saturated 
deposition (see Section 5«2o2) the spin torque is about 1 
dyne-cm which Is quite apprecialbe for long term spin 
control because the life span of this thrust is about two to 
three half “lives of the alpha-emi ttero 

6olo5 Erecting rate 

During the despinning period/ the erecting shell will 
nearly attain the equilibrium temperature distribution/ and 
the nutation will be damped out by the annular viscous 
damper (see Section 2o3ol)o Assuming the transient 
precession and erecting torques do not affect the attitude 
significantly/ the Initial motion of the spacecraft is 
simply a slow retrograde precession about the sun axiso 
Since the precession torque and the spin remain relatively 
constant/ the effect of the erecting torque on the 
spacecraft can be found from the M 2 expression of Equation 
(6o3)o The resultant erecting rate In terms of the erecting 
torque is 
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©•“Mj/ C Ci^ “2 It \|/cos@) 

For a practical design^ the net precession torque (algebraic 
sum of the solar pressure torque and the parallel torque 
component from the erecting shell) is very small ^ hence the 
precession rate Is usually very much smaller than the spin 
rate« The erecting rate Is further simplified to 

e— Mj/Cl^oj^) CS.10) 

A shell ( thi ckness<*0o01 cm^ radlus»25 cm^ 1ength*25 cm) will 
yield an erecting force of OoOl dyne at the optimum spin of 
O 0 O 37 rad/$ec^ so an erecting torque of 0o25 dyne-cm is 

achlevable« As a result^ the erecting rate for the bare 

■"6 

shell is about l.»3*10 rad/sec and this rate will decrease 
as the spacecraft precession cone angle is decreased because 
of the reducing solar input to the shell surfaceo 
Therefore^ It will take on the order of two weeks to 
stabilize the spacecrafto 

The main drawback of the above scheme Is the 
vulnerability to small perturbation torques because of its 
low spin rateo A radi o I sotope°° heated shell can resolve this 
difficulty^ For a surface heat flux of 0»2 watt/cm^ the 
optimum spin rate Is 0ol4 rad/sec and the corresponding 


erecting 

torque i s 

about 

0o27 dyne-croo 

Due to 

this 

large 

Jump on 

the spin. 

the 

erecting rate 

i s down 

to 

4ol*10^ 

ra d /sec 

0 Consequently, 

the stabilization time 

increase to 


about one month o 
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For short-terii! stabll I zatlon^ the subUmlng raaterlal 
coated shell Is really the best performero Even at a spin 
of rad/sec (way above the optimum spin), the erecting 
force is nearly 20 dynes and the erecting torque Is about 
^^50 dyne“cmo The erecting rate Is 3*10 rad/seG which can 
damp out the precession In a few hours at mosto 



6o2«l Introduction 

In general, any of the above erecting shells can be 
made to work for a solar probe, like the present Sunblazer 
satellite. If the spin Is controlled actively., For 
Instance, a stepping motor may be used in conjunction with 
an attitude sensor to switch the spin torque by turning the 
microthrusters such that the spin of the spacecraft can be 
Increased, kept constant and decreasedo 

Disregarding the present Sunblazer att itude-control 
design philosophy, the erecting shell and the microthruster 
are able to fulfill the stabi 1 Izatf on requirements for 
future solar probe missionso Three conceptual llght“welght 
attitude control systems are described in the following 
sect Ions o 

6*262 Passive systems for inferior orbit probes 

For Inferior orbit probes, the solar Input is quite 
adequate so the pure reradiatlve torques from a thin 
metallic shell and a restoring vane (reradiatlve Internal 
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surface) can be acceptable in the following passive system 
designo 

The chopped°°cone micro thruster is designed to have 
enough subliming material such that the spacecraft is very 
rapidly despun shortly after launGh to the vicinity of 
optimum spina in this low spin range^ the erecting torques 
from the thermal reradiations (summation of the effects of 
the shell and vane) are most effective in reducing the 
precession cone angle and even at zero spin^ the reradiatlve 
force is still capable of damping out the resultant 
libratlonal motion (2)o During this erecting period^ the 
small recoil thrust due to the radioisotopic emission of 
alpha particles from the depleted microthruster despins the 
spacecraft to zero and then starts to spin up in the 
negative direction so the effects of the perturbing torques 
due to the interspace magnetic fields meteorites 
bombardardments and outgassing from leak-joints can 
gradually be reducedo This is by no means a fool-proof 
system because the exact dynamical behavior of the 
spacecraft near zero spin is still not well understood, 

6,2,3 Semi-passive system for Inferior orbit probes 

The disadvantages of the passive system are mainly the 
slow response and the inflexibility. The above passive 
system can be improved by installing a motor which^ on 
command from the attitude sensor^ is able to switch the spin 
moment from positive to zero or negative and vice versa. 
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One schefiie which can make good use of this 'switching' 
motor and attitude sensor combination is the foliowingo A 
cylindrical strlp-gap vane as shown in Figure 6,4 with the 
internal surface coated with subliming material Is placed at 
the rear ©f the spaceerafto At any instant^ the attitude 
sensor detects a deviation from a preset precession cone 
angle^ the microthrusters will be positioned properly by the 
switching motor to despin the spacecraft toward a 
predetermined spin rate where the undeslred positive 
precession torque due to the parallel component is 
insignificant yet the erecting torque due to the normal 
ccKQPonent Is still quite appreciable to stabilize the 
spacecrafto As soon as the desired attitude is attained 
again^ the microthrusters will be switched so the spin will 
go back to the original steady state value at which small 
distrublng torques are not likely to upset the attitude of 
the spacecrafto in the event the attitude is thrown off^ 
the above procedure is repeatedo Since subliming materials 
on the vane and In the micro thruster must be present 
throughout the mission^ the sublimation rate during inactive 
period should be less than one cm/year « As a result^ the 
steady state surface tetiperatures have to be less than 225" K 
In order to reduce the weight penalties for the subliming 
roaterialso The sublimation rate and pressure as functions 
of temperature are shown in Figures 6o5 and 606 for 
naphthalene and biphenyl 0 For long missions^ biphenyl Is a 
suitable material because of its low subliming rateo 
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Intuftively^ the thrust level Is dependent on the surface 
radiation properties as well as the gap°°strlp ratlOo The 
cylindrical vane shown in Figure 6»4 should reduce the 
steady state fuel loss through continuous subl imationo 

6o2o4 Active systera for superior orbit probes 

For probes designed to explore the electron density 
and Interspace laagnetfc field beyond the Earth^s orbits and 
up to 50 AU distance such as pioneer scouting satellites for 
the Grand Tour Riission^ almost all present systems have to 
be raodifiedo Because of the 1/r^ attenuation on the solar 
flux^ the electric power generation will have to turn to 
nuclear devices such as the radioisotope thermoelectric 
generator CRTG) or the radioisotope thermionic devIceSo 
Also^ as result of long transmission path^ a high”galn 
antenna Is necessary which makes the attitude control 
requirements quite formidableo 

Although^ the erecting shell and the restarting vane 
will lose their appeals^ the microthruster can still be used 
to supplant their original functions (generation of the 
erecting torque and the negative precession torque and the 
despin torque) o This task Is accomplished through an active 
coupling of the mlcrothruster and the attitude sensor 
systemo 

The sequence of Interactions to stabilize the attitude 
Is the followingo From Figure 6o7<? the spin rate and the 
relative angular positions of each pair of microthrusters 
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with respect te the precession axis (loeo^ the angle between 
the line connecting each pair of microthrusters and the 
precession plane or the plane of the sun^polnt axis and the 
body axis) can be determined by the attitude sensors « Then 
a trigger signal delayed by a preset spin angle starts the 
turning of each pair of microthrusters at the same rate as 
the spfPo in each spin cycle all three moments are 
generated as shown In Figure 6o7o The direction of turning 
will only affect the erecting torque and the precession 
torque but not the spin torque which is always averaged to 
zeroo The time~de1ay on triggering the turning will affect 
the waveforms of the precession torque and the erecting 
torque as shown In Figures 60 S and 6o9o 

The maximum averaged precession or erecting torques are 
only a half of the maximum spin torqueo These torques are 
independent of the spacecraft attitude^ in particular the 
precession cone angle^ with an exception that as the cone 
angle reduces to zero the attitude sensor will not be able 
to trigger the turning properly (assuming the sensor Is 
attached at the periphery of the spacecraft) o Since the 
precession is a function of time ^ the attitude sensors may 
have to trigger frequently for Instance^ one trigger for 
every two or three spin cycles during the stabilization 
pertodo Once the desired attitude Is at ta I ned.^ all torques 
can be reduced to zero by turning the microthrusters such 
that the spin torque Is zerOo Whenever the attitude sensors 
detects large deviations^ this stabilization system will be 
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activated again<, 

In principle^ heat°’pipes can be used to control the 
heat flux supply for sublimation by mechanically causing the 
heat*’plpes to contact a heat source whenever heat is 
requiredo in so doings the fuel requirement can 

significantly be reducedo Alternatively^ the long-lived 
radioisotopic recoil momentum can be used for generation of 
thrust^ however^ the stabilization time will have to 
increase because of the low thrusto The stabilization time 
can be shortened if the spacecraft is despun whenever 
stabilization is called for^ and as soon as the 
stabilization is attained^ the spacecraft is spun up to a 
steady state spln^ at which the spacecraft is relatively 
immune to external torques o 
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FIG 6.2 DEFINITION OF EULER'S ANGLES 

lo 

1 1 ’l i I 



FIG 6.3 GRAPHICAL REPRESENTATION 

OF Eq(6.4) 
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FIG 6.4 CROSS-SECTIONAL VIEW 
OF STRIP GAP VANE 
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PRESSURE dyne / cm* 



6 





SUBLIMATION RATE cm /sec 


FIG 6.6 SUBLIMING RATE 
vs TEMPERATURE 
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FIG 6.7 

RESULTANT THRUSTS FOR SYNCHRONIZED 
TURNING OF MICROTHRUSTER (trigger 0 = 0 ) 
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FIG 6.8 

RESULTANT THRUSTS FOR SYNCHRONIZED 
TURNING OF MICROTHRUSTER (trigger 0=9Or) 



variation of 


thrust In a spin cycle 
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FIG 6,9 

RESULTANT THRUSTS FOR SYNCHRONIZED 
TURNING OF MICROTHRUSTER (trigger 0=135°) 

so[lai| Injput 



Variation of 




CHAPTER 7 


CONCLUSIONS AND RECOMMENDATIONS 



The eoncluslons reached In studying the appllcatlens of 
radioisotopes and subliming materials for attitude control 
of small solar probes are listed belowo 

(1) A1pha°°emi tting radioisotopes with moderate 
half-lives are much superior to beta and gamma emitters in 
applications for attitude control » These isotopes can 
readily Increase the optimum spin of an erecting shell so 
that the spacecraft is more Immune to disturbing torques., 
Furthermore^ they are Indispensable in supplying uniform 
thermal energy to a subliming microthruster for fast spin 
control 0 Additionally^ they are useful to generate small 
recoil thrust for long term^ fine spin control o 

(2) Organic sub liming materials are found to be very 
suitable for use in magnifying the erecting torque from an 
coated erecting shell by several orders of magnitude over a 
wide range of spin In comparison with a bare thin metallic 
shello As a result^ the erecting shell can still function 
effectively even at a relative high spin rateo 

(3) A miGrothruster fueled with subl iming material can 
convert almost all the thermal energy imparted to the 
sublimed molecules Into useful thrust by proper design of 
the nozzle cone angle and lengths 
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ik) At a given sublimation rate^ the thrust can be 
further fncreased if a portion of the nozzle skirt Is heated 
radio I so top I ca 1 1 y « 

C5) The applications of ffllcrothrusters and/or erecting 
shell are found to be attractive for stabilizing the 
Sunblazer probeo With proper feed°°back control mechanisms^ 
the appl leatlons of the above devices are able to be 
extended to a wide variety of satelUteSo In particular^ a 
pair of microthrusters can produce torques in all three body 
axes to control the attitude of any light spacecrafto 

Investigations on other required studies led to the 
following conclusionSo 

<11 In many cases of Interest^ the transient thermal 
analyses of erecting shells and microthrusters^ can be 
carried out by considering only zero or one spatial 
dependence without appreciable loss of accuracyo Numerical 
instability Is found to limit the time step for ADI solution 
of a subl imlng ralcrothruster where both radial and axial 
beat conduct Ion effects are accounted foro 

(2) The transient period for an erecting shell Is short 

(less than an hour) if the attitude of the spacecraft is 
fixed o However^ during stabilization^ the attitude 

constantly changes^ so the transient effect should be 
included In detailed dynamic analyseso 

(3) The numerical technique (modified Monte Carlo) 
developed for the non'^obst ruction factor calculations of 
molecular flow is very flexible^ accurate^ and can Include 
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complex physical phenomenao The results were applied to 
study the evaporation process and found to be usef u1 o 

(4) The above numerfeal technique was adapted to 
analyze the gas dynamic aspect of the subliming 
mlcrothrustero The result clearly paved the way to an 
Improved design of the ralcrothrustero 

7^^ Recommendations for further studies 

The present study is by no means completeo Several 
areas are ©pen for immediate further research o 

(1) The method of calculation of the molecular 
collision effect on the non“obstruct ion factor needs 
refinement^ since the evaporation coefficients derived from 
existing experimental data of mass flow experiments appear 
to be low thus suggesting that the molecular collision 
effects used in obtaining the evaporation coefficients have 
been underest I matedo 

(2) The technique developed for non-obstruct i on factor 
calculation should be applicable to shielding calculations 
involving neutron and gamma ray streaming through a gap or a 
beam-port t, 

C3) The assumption that the sublimation surface 
reflection coefficient is one minus the evaporation 
coefficient can be tested experimentally as outlined in 
Section 3o2o3o If it is verified^ the evaporation 
coefficient theory will be placed on a firmer basOo 
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(4> In the Gh©pped«eone micrethruster analysis^ the 
sublimation surface Is assumed flat^ but in reality the 
surface becomes coneave« The effeGts of this curved 
subliming surface on the thrust are yot to be investigated, 
C5) Utilization of molecular dissociation may offer a 
means for a marked Increase in mlGrothruster performanGe/ 
but a radical change In the design proposed here would be 
required in order to provide the few electron-volts of 
energy required per molecule to induce molecular 
dissociat ion, 

(6) Detailed analyses should be carried for the dynamic 
behavior of solar probes which Incorporate the various 
erecting shell and microthruster designs proposed here in 
order to determine the specific design requirements for 
specific missions. 
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APPENDIX A 


APPROXIMATE ANALYTIC SOLUTION OF t-6 
PARABOLIC PARTIAL DIFFERENTIAL EQUATION 

The partial differential equation Is of the form 

Aff = tCT^+DFCe.t) (fl.j) 

where symbols are defined In Appendix is used for Thin 
Shell Approximation In Section 4<,3o4o At steady state^ the 
temperature distribution has the approximate form of 

“T( 6, t) - [l + M 

By trigonometric expansion^ T (S^t) becomes 

TCe.t) = [ 2 + M, coscMt-t ej+ M^siNcuJt+e)} 

and 

I^X = 7 ;^ cOSCi^t-tO) -Mz S/A/CaJt-t9)] 

H- 4-M>C0SC^t-ta) + ‘^MxSlNCi^M)lC^-^) 

The solar Input F(B>t) Is a half-rectlf led sinusoidal 
and can be expanded as 

fCB,t)= CA.7) 

Substituting the above quantities Into the governing 

differential equation and neglecting the higher order termso 

+ MzCoscu)t±6)] = 

-BT:^j[M,C0SCu}tt8)fMs SINCwt-i-6)] -t- + 

4-Mi cos Co)t-t &)+4 M2.s//^coJt+0 )] -f X>[-^ +^SWCM*ej] 

CM ) 


CM) 

CA3) 

CM) 

Cf\-S) 


185 



Since the above identity Is forced to hold for all t and 9, 
the coefficients of the slne^ cosine and constant terms must 


be zero^ toe«^ 

C + 3>/ir -o ^ 

3 Ta^^M, + =d rRJl) 


The solution is 




So that 


o-5fi Ta^^ iX) j> 

0.5 I>(-BTav^ t4cT.KVg^) 

- B u)"- IGc^t/v^ 1- BBC 


M.» 


M M/ + m1 


( A. IH) 


(A. 13) 


()R. /f) 


. Cfl. /5) 


lo Cos 4> 

“ 2^r if>^c‘oy^i- KVr^ +s-^Xv%+ie.^X.v^]^" 

< - M. X 

i = tan (— ) 

= tan-( 

For subliming material coated shelly a heat loss term ET 


EXP(Ae”E^T) Is added to the differential equationo The 
expansion for this term is 

r p ( ft - Bc/t) - Tcvi'g) t'^I^ ^ 

fi i-( -M£/2) stn Ct^t-t©) 

+ f -H,/2 ) cos Cu)t -»■ e)J (R-l"]) 
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Substituting back^ ^ M, and are solutions of the 

following equations 

c + ^/n 

PI (/0M2 = ' BTav^M, t £ * 

exp(R^- ) 

t- (AvR) 

= 'BT^>'gK)2 +4 CMi T1^-^ + 

£ ^ ) ( BcM zIx^^/^ - Mi /2 ) 

The results for M and are given in Section 4«3o4o 
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APPENDIX B 

APPROXIMATE ANALYTIC SOLUTION OF OPTIMUM SPIN RATE 
In Section 4<,3o4^ the apprexiroate analytic optiimmi spin Is 
obtained as following. At the steady state^ the teBiperature 
distribution Is approxlmatel y given by 

= (B. i) 

In an Inertial fIXed coordinate^ the teroperature 
distribution Is fixed as shown in Figure B.! and Is reduced 
to 

T ( ) [ 1 t M S i h ( 9 t ( g 2 ) 

The net force component In 9 direction (or normal to 
solar ray) is proportlondl to the following integral 

.2TT 4 

co>s eda (B. 3) 

^ 0 

Taking the derlvatlye with respect to and setting the 

result eoual to zero^ the result Is 

4TT2 w‘+f SlTj vt TCj + airf TTj 

+ (6-irf iT 2 - 4WiTf, ]<*)*- 3 n-j=o ( B- +) 

where parameters are defined In Section 4.3 .4 and Appendix 
Mo This is a 6th order equation from which a positive real 
root Is to be solved. There should be only one real root 
which Is the optimum spin speed for a maximion erecting 
force. 

The method of ^nalytlc solution to this 6th order 
algebraic equation follows. Transform the equation to 
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2 

eubicCby setting to=y^ then eliminate the quadratic term by a 
linear transformation (y=x+K)^ where K“~(2iTjTri7^+5Tr,V2 
The cubic has the form of 

PX + CL = O (B, s ) 

Now^ the three roots are by Cardanos Formula ijjJJ 
XI ^ ^ +B 
KZ - cOR + 

)(3 = (B.6) 

o^\\Q.re. 

«^=L-1 '•-/3 OA 

fl’ = - Q/e +Ja74 1-P%7 

-a/^ -JqV 4- + P%7 fB- 7) 

Since the discriminant Is negative^ all three roots are 
realo From the DesCartes® Rule of Sign (H)# the original 
6th order polynomial and its transformed cubic equation have 
one and only one variation of slgn^ so that the equation can 
have at most one posIt|v|e real rooto On the other hand^ the 
cubic equation can have at most two negative real roots^ 
hence the original equation can only have one positive real 
root which Is the desired optimum spln« The closed form of 
this root given In Section 4.^3 » 4 was obtained by means of 
the complex variable identities manipulation of Equation 
(Bo6)o 

The numerical results from the analytic expression 
Equation (Bo 6) have been checked against the direct 
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numerleal solution of the original polynomiaVo The six 
roots Indeed contained only one positive root which agreed 
to four or five significant figures with the closed form 
solution., 
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FIG B-1 STEADY STATE 
TEMPERATURE PROFILE 
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APPENDIX C 


GRAPHICAL EXPLANATIONS FOR NSHELL SOLUTION 


For a given configuration of the erecting shelly the 
steady state temperature profile around the shell with 
respeet to the sun is Independent of ttme^ however^ the 
transient behavior Is very much dependent upon the Initial 
conditlonse The steady state temperature profile is 
characterized by the maximum T^^, the minimum T^,„ and the 
phase lags f and ^ as shown In Figure Colo To a good 
approximation^ their relations are 

= cc.l) 


1 . ^ 


C c, 


Case 1 <T 




I. 


i I I 



where t^ Is the time required to rotate to Lj. and 'TV is the 
time required to rotate to L^ from the positions shown in 
the figureso 
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I f, 

However^ If is very mueh larger than ©sell 1 at I ©n 


will disappear^ instead the transient distribution will look 
1 1ke a slanted stairo 






APPENDIX D 


RECOIL PRESSURE OF PHOTONS# NUCLEAR 
DECAY AND SUBLIMING PARTICLES 

(O Solar pressure on a projected surface is 

Psu.n ^direct solar pressure + reflected solar 

pressure 

) I, /c 

- Iglc. * il * C t~ & ) 

where 1° is the solar flux at a distance 1 AU 
G is the speed of light 
& is the emissivity 
( i i ) Thermal rad i at i on pressure i 

For a diffuse surface the total emissive power is E^irlo 
The time rate of change of normal momentum is 

* 2 . ^ cj^d ^ j ^ ^ ^ ^ 

Vi 

pj, = 2 E/c ^ 5 4* ^ 

^ o 

= 

(I.. 3) 

Cilt) Nuclear disintegration pressurei 

The recoil mcMnentum due to a relativistic partiGle is 
E/C*(l+2mC^ / e/\ The number of disintegration per unit area 
is 
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(p.^) 


Nsdislntegrat !on per Watt * Watt per area 
=3o7*lo'“ * 1, 686*10* /E(tnev)*WPA 
The dlstntegratien pressure is then 

0,832. *10^ ^pR * ^ / E iV- S ) 

*This expression is good only If the deposited layer 
thiekness Is much smaller than the alpha range In the 
material/ i«eo/ self absorption is negllglb1e» 

Civ) Sublimation pressurei 

The recoil pressure due to sublimation Is quite 
complex/ but as a definition (J.) 

Ps ' ^ '^'^1 Z. 

~ 6 . (P.6) 

where Py is the sublimation vapor pressure 
o( is evaporation coefficient 
J is the subliming temperature 
* The expression is good only if the reflection is specular 
and con Stan to 
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APPENDIX E 


DERIVATION OF THE RELATIONS BETWEEN 
RECOIL PRESSURE AND REFLECTION COEFFICIENT 

Let the interactions of vapor moVecules with solid 
surfaces^ i«eo^ the reflection (dif fuse+specular) 
coefficient as a function of the Incident angle speed v 
and surface temperature T be r(S#v^T)o The fraction of the 
molecules which do not reflect are assumed to be condensed. 
The fraction which does not reflect Is assumed to be 
condensed. At kinetic equilibrium^ the fraction condensed 
is equal the fraction evaporated. Therefore^ the total 
number of molecules evaporated from a plain surface is 


proportional to 



r c©,V.T)Jf (6 / O dv 9 


where fC0^v) is the number of particles reaching a unit 
surface area per unit time and solid angle 
with speed between v and v+dv 
At thermal equilibrium^ the vapor*‘pressure acting on a 
surface consists of the contributions from evaporation^ 
condensation^ diffuse and specular reflections. That is 
Py- Pe+ Pc Fd t Ps 
Pe - Po 

Pet _ fye^e B'vcose 
Pc ^ ® ■" r V cose 


rE - 1 ) 
CE- 2) 

(£.3) 


P5 V CoS ® t CoS B 

p6 V* Go's © rv cos© 


(p4) 
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whftre the barred quantitfes are the Maxwellian distribution 
weighted averageso The reflection coefficients are 

functions of the surface temperature^ the Incident angle and 
speedo The total reflection K' I s the sum of the specuiar s 
and the diffuse reflection 5 ° Substituting the above 

expression into Equation (Eol)^ the relationship between 
and Pe is 


n 

Fb 


e + 


rv<^o 5 a ir S 'V(u>^B -i- s\/c^ss 


(B.S) 

vc-flse-t'V'Cos5 '' 

Only if 8 is zero and a is independent of both the 
incident angle and speed will the above expression reduce to 

Vc/-r) 

- s ) - oi C B. S) 

where c< is the evaporation coefficient 
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APPENDIX F 


GEOMETRIC RELATIONS FOR SPECULAR 
PARTICLE REFLECTION IN SPACE 

£aJLjExJ,,lfl.drJ.s,a,l MaT 1 

A particle ©riginating from P, (x, ^y, ^z, ) and Impinging 
the wall at P^ (x„^yo^Zo) will be reflected and will pass P^ 
ix,y,z)o The other point lying on the normal of the surface 
of reflection Is P^^CO^O^z ), Since are In a 

same plane and Pj Is a mirror image of P| # a plane equation 
and two sphere equations can be written to solve for P^ 
(x^y#z)<> The results are 

X~ iJ/r" 

^ = <F. E) 


z= 2 Z, ' 2^ 


(f.3) 



on the wall^ the Intersection point has to be solved from 
the line (pass! 
simul taneousl y, 


the line (passing through P^^l^) and the cylindrical surface 


£.aij£h,eftefidj;.,gns. , MaJ.]. 

The specular reflection from a chopped cone wall Is 
more complex because the surface normal is not readily 
obtainableo The procedure to calculate the reflected beam 
Is similar to that of a cylinder^ but the resultant 
expressions are very comp1leated« They are used directly In 
Subroutine CONREF In prograra CONNOF (see Appep^lx Ko2)o 


198 



APPENDIX G 

DERIVATION OF PARTICLE FLUX DUE TO SUBLIMATION 
The angular mass flow rate from a differential area dA 
for molecules having velocity In v and v+dv Is /4 tt ) vdn^ 

cos&dAo The number of these molecules streaming through a 
unit sphere is C<=^/4TT)vdn^eos0CdS/ ,C^)dA per unit tlmco The 
average residence time of these molecules {nsld^ the unit 
sphere is 4r/3Vo At any Instant^ the number Of these 
molecules stationed In the sphere Is ( ®^/4ir )vdnu- c©s8dS/je^ 
)(4r/3v}dAo The flux or track length due to these molecules 
Is (Q(r^/31^)vdnvrCos6dA or (‘=*/4 tt 1’^ ) vdn^ conBdAo The flux 
from molecules of all speed is just (p^n/(4TTl*') )vcps6<iA Or (ok 

— 3 

nvz/4ir I )dAo Summation of contributions Of all the 
elemental area dA will yield the flux at the location of the 
sphere due to the sublimation onlyo The corresponding 
particle density is n»^(dA^nz/l^ )<. Contri butlons fr^ 
channel wall reflection and reemission are not Included In 
the above summationo 



APPENDIX H 


DESCRIPTIONS OF COMPUTER PROGRAMS FOR 
THERMAL ANALYSIS OF THE ERECTING SHELL 

This prograra (STRIP) solves for the temperature 
distribution around a shell by neglecting the effects of 
thermal conductlonso The shell divided into UM strips., At 
each time step (Is fixed by JM)^ the temperature Is set to 
be uniform; however it can be very easily modified for any 
initial distributionp The maximum^ minimum temperatures and 
phase lag at each time step are computed and reGorded for 
plotting at the end of the computation determined by TMAX in 
seconds p 

Subroutlness ARRAY^ PLOT^ FORCOM^ CURVEp 
Functions! FUN, DIFQ2, FH, FN 
Compilers WATFOR 
Execution times Op 7 minute 
- She 1 ,1 App rox I wa t jjan 

The complexity of the coding Increases significantly by 
including the tangential heat conduction effectSp The 
coefficients of the difference equation for each strip are 
computed and stored in the coefficient matrlxp The number 
of strips, UM, is preferred to be odd and less than 50, For 
each time step, the matrix is solved by Gaussian reduction 
method, the reradlative force components, maximum, minimum 
temperatures and phase lag are computed and recordedp The 
Interaction is terminated as soon as the NMAX (<200} time 
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steps or the steady state has been reachedo The stored 
informations are plotted out and the steady state 

temperature distribution is expanded numerlGaily into 
Fourier sine and cosine serleSo The approximate analytic 
results are computed for eomparisono 

Subroutiness RECOIL^ FORIT^ RGRAPH^ PLOT, 

ARRAY, CURVE, APPROX, SIMQ6 
Functions SiNHR 
Compilers WATFOR, FORTRAN G 
Execution times lo5 minute 
H.3 Thick Shell Approx imat_LQn 

By considering the radial heat conduction alone, the 
partial differential equation remains one dimensional o 
However, for each time step, the equation has to be solved 
as many times as JM (number of strips) <> The difference 
equation for each concentric layer yield a tri diagonal 
matrix which are stored in AA and solved by the fast matrix 
factorization raethodo The surface temperatures are plotted 
out occasionally both in rectilinear form and in polar form<, 
The maximum, minimum, and phase lag for the internal and the 
external surface are plotted out every 100 time steps » 
Subrout iness 0UTPT2, STRIP, MATFAC, PHASE, MAXMIN, 
ARRAY, PLOT, CURVE, 0IFQ2, POLAR 
Functions s FUN, SINHR 
Compilers WATFOR 

Executfon times 3 raino/100 time steps 



H.4 Optimum Spin Rate 

This program Is relatively sitnpleo The coefficients In 
the partial differential equation are computed firs to The 
average temperature Is solved by systeni subroutine RTWI 
(using the Welgsteln^s method) <> Then the coefficients for 
the sixth order algebraic equation are computeda The six 
roots are solved by another system subroutine PGLRTp The 
numerical value of the analytic optimum spin expression is 
also ealculatedo 

Subroutiness RTWI<? POLRT 
Functions g FCT^ FCN 
Compiler; FORTRAN G and H 
Execution time; 0o3 minute 



APPENDIX I 


SPACE CURVE TRACED BY CONIC-PARTICLE 
BEAMS ON CYLINDRICAL TUBE 

I »1 Inlet Sublltaatlon 

The equation for a conic surface Is obtained by 
rotating a straight line about an axlse This conic surface 
can describe the sublimation from a tube in1et» 

A line In X-Z plane is 
Z=x*tan6'-h*tan8 

The surface generated about x=h is 
Z “ tan^ S C ( X- h f + y ) 

The equation for a cyl indrlcal tube Is 

•2. % Z 

X +y =R 

With the above relations^ the parametric equations for the 
space curve of intersection are 
x«(R^-t1/2h 
y«(l^ -x^ 

z “i/a 

Z=(h + 1) ' tan@ 

where t Is the parameter 

These three equations yield the coordinates of a point on 
the space curve ©f intersectlono 
l„2 Channel Wall Reemisslon 

in this case^ the reeroission Is originated from the 
tube wallo Ti>e conic surface 
(R-x)^ »tan^@( (z-z, f +y^) 
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eoiaes frora the revolution of the line 
R“ 2 «<^-z^ )tan9 

Where (R^0^z») is the apex of the eone and 8 is the 
angle between the generating line and the z axiso The 
oylindrieal surface of radius R is 
X +y “R 

The resultant parametrie equations for the space curve of 
intersection are 

Z= Z„ 

(R~x)^ ®t*tane"(x”R^ ) tan^e 
s. ./2. 

y«(R “X > 

Again/ t is the parametero This result is useful for 
computation of the fate-probability of reemitted parti clesa 
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APPENDIX J 


DESCRIPTIONS OF COMPUTER PROGRAMS FOR 
THERMAL ANALYSIS OF THE MICROTHRUSTER 

J>1 Cvl Indr leal Tube—Exol id t Method 

Since the temperature on any mesh point at any time 

step depends only on the temperatures of previous time step^ 

The transient solution is obtained by marching-out In timoo 

However^ due to sublimation^ the extreme mesh always 

decreases in sizeo As soon as the numerical stability 

criterion is violated^ the mesh net is rearrangedo Toward 

the end of sublimation^ the rearrangement will introduce a 

large dlsturbanceo Aise^ the time step has to reduce 

frequeniy to insure stabilityo The surface and the bottom 

temperature and the fuel thickness are displayed graphically 

by PLOTo 

Subrout I ness OUTPTl^ CURVEl^ ARRAY, PLOT 
Functions G 
Compilers WAT FOR 
Execution times 2^5 mlno 
d.2 Cvl tndrieal Tube— I mpl I el t Method 

Unlike the explicit method, the temperature 
distribution in the fuel at any time has to be solved 
simultaneous lye Since the Goefficient matrix is 

trldlagonal, a fast matrix factorization method is adoptedo 
The initial temperature distribution is arbitrary set to be 
linear 0 The method is numerically stable for any mesh sizoc 
However, toward the-end, the mesh net has to be modified to 
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twe**p©]nt differetieeo A t«mporary disturbance did show up 
In the surface teRiperatureo As usual ^ the results are 
plotted oUto 

Subroutines? MACFAC^ OUTPUT^ CURVEl^ PLOT^ ARRAY 

Comp Hers WATFOR 

Execution times 2mlno 
J . 3 Chopped Cone- ADI. J<et hod 

In the analysis of a subliming microthruster^ due to 
the cylindrical symmetry^ the temperature In a chopped cone 
ralcrothruster is a function ©f r and Zo The physical 
problem Is best described by the following diagram^ Figure 
Jolo The presence of sublimation^ which Involves both heat 
and mass transfer# prevents the temperatures from reaching 
steady states To apply the ADI technicue# the Gone cylinder 
is divided Into convenient mesh as shown in the 
cross-sectional diagramo Difference equations implicit in 
r-direction are required for each row as well as similar 
difference equations implicit In z-direction for each 
column^ The diagram shows that there are only 13 distinct 
difference equations for each direction (ioeo# temperature 
at other me;sh point can be expressed by one of ti^e 13 
equations)o Their formulations are strictly based on energy 
balance (first law of thermodynamics) and with proper care 
no problem should arise for difference equations In the 
center reglono 

The general heat conduction Is Q«-K*S*D* To For 
rectangular mesh S*D Is OoSAAx (note? D Is depth# S»0o5*cot& 
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y where B Is the angle subtended by the heat flow path 

The temperature at any mesh point can be described by 
one of the 26 (2*13) difference equationso At each time 
step^ the temperature distribution In the plane Is obtained 
by solving the temperature in the line repeatedly until the 
plane Is eoveredc. The direction of the lines alternate for 
each time stepc The rep resen tat i ye temperature profile is 
shewn In Figure do2o The iteration in z-dl recti on Is 
performed In Subroutine ZITER and the Iteration in 
R-di recti on is performed In Subroutine RITERo 

Sub routines s RITER^ ZITER^ MATFAC^ PARAS, FARAD, 

RADFOR, CURVE 
Functions SU6M 
Compilers WATFOR 
Execution times 1/3 real time 



FIG J.1 SCHEMATIC OF CHOPPED 
CONE MICROTHRUSTER 




45 degree cone; radfus^OoS cm; time step»l second 

fuel»naphtiialene; thicknesf »0o95 cm; emlsslvlty«!i,0; heat source=0o75E 07 erg/cc-$ec 
cap mater I al»6e; thickness <^0 0 05 cm^ emlssivity*0o5; heat flux= 0 o 5 E 07 erg/cc-seco 
Ihitidi temperat4re®273£,0 fc (uni form) 
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APPENDIX K 


DESCRIPTION OF COMPUTER PROGRAMS FOR 
RAREFIED GAS NON^OBSTRUCTI ON FACTOR 

IL I .Cxl I ndrJji.i.1 iMfeg,.. torafes ,t r»s.t Easlgi: 

After Initialization of the necessary parameters^ the 
relative amplitudes for all angular groups in the spherical 
cosine are computed and stored^ then a probability ’table’ 
for reemi tted molecules Is constructedo Now^ molecules 
leaving the point sources are traced success! ve1y« 
Symmetric properties are utilized In order to save execution 
time whenever possfbleo The molecular collision effects are 
estimated by calculating the molecular flux and density in 
the channel first and then approximated the fate of the 
collided moleeuleso Finally^ the fate of the adsorbed 
molecules are computed with the aid of the probability 
’toble^o The total exit molecules to the initial inlet 
molecules rati© Is the desired mass non-obstruetion faetor^ 


Subrout i ness ANGP^ 

COLISN^ 

XY, 

MOLCOL, 

DELAY, 

RE EM IT 

XYZ3<, 

ITER, TABLE, 

RGRAPH, 

XYZ2, 

XATER, 

RFLUX^ 

DFLUX, 

REFLET, BISRCH, SCATTER, 

CONE 







C©R«)llers WAT FOR 
Execution times 1 ml no 

The basic techniciue In this program is essentially the 
same as the previous ©neo However^ the molecular collision 
option Is replaced by the subliming surface reflection and 
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the thrust coefflelento Also^ the mesh e©ne Is divided Int® 
tna®y equal azimuthal segments which are different In the 
cylindrical case where the mesh cone is divided according to 
the Intersection of the mesh cone with the wall ring segment 
rlngs« For 90* cone angle (loe*^ cylinder)^ the result .can 
check against those of previous calculat lonso 

Suhroutlness HIT^ CONREf^ CONCON^ ANGPR^ CONCEN^ 

BISRCH^ CONEMI 
Functions? CNITER^ CENTER 
Compilers? FORTRAN H and WATFOR 
Execution time? 0o7 mino 
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APPENDIX L 


SUBLIMATION PRESSURE 

In caleulatlng the subllfflatlen heat flux or mass flew^ 
the pressure term always ariseso The pressure depencence on 
temperature Is derived belowo The free energy ehanges is 
constant for both the solid and the gas phases at 
sublimation for small changes In P and T^ therefore 
dG»Vv.dP-Si,dT-Vs dP»Ss dT 
or dP/dT»^S/^V 

M 

if the latent heat ©f sublimation is ^ ^S is simply 
Assuming the Ideal gas law, the following approximation 
can be obtainedo 

dP/dT»Hs/C T^V) 

“Hs/TVg 

dP/{PdT)-Hs/(RT) 
dlnP/dT«H5/(RT' ) 
lnP»“(H5/RT))+C 

. A -B/-r 

P*exp ( + c ; “10 

This is the common pressure*- temperature relationship for 
sublimation^ A and B are empirical constantSo 
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APPENDIX M 


NOMENCLATURE 

A ,B -*“sub1 linat i©n eonstants 

AU —“astronoinleal unit (radius ©f earth ©rbit) 

A P© 

B K/R* 

e —speed ©f light; heat eapaci ty 

C — (-6-^/^r) 

d -‘-diameter 

D — -ioeos4/'“r 

E — (-H^Ft/ClJolfe^r)) 

G — RlHs/^r; mass flew rate 

H — angul ar momentum 

Hs —“heat of sublimation 

To —moment of inertia 

Ivp ' — -solar eon Stan t (heat flux) 

I >spee i f i G i mpul se 

j; —“thermal conductivity; BoltzroannTs constant 

K —non-obstruction factor 

1 —channel length 

L —subliming fuel height 

ID —particle mass 

M — -molecular weight; temperature amplitude factor 

— -moment or toraue component In I -direct I on 
— particle density; neutron 
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p 


—pressure; mementtBii 


q ---heat source 

q" q -o-heat flux 

r^ - — alpha range 

R — eyUndrlcal radius 

RIm* —radioisotopic heat flux 

s —specular reflection coefficient 

t — tirae 

T —temperature; kinetic energy 

— -average temperature 
V - — particle speed 

W* — molecular weight 

oi, — ^evaporation coef fee lent; alpha particle; absorptivity 
)r —reflection coefficient; spectflc heat ratio 

5 - — diffuse reflection coefficient 

fr °°eml ssl v! ty 

p —density 

a-' — -Stef an-Bol tzmann constant; moleeular cross-section 

« 

4' —precession 

Np —maximum temperature phase lag 

f —angle between figure axis and the sun- line 
co ...gplq rate; angular speed of body 

^ —minimum temperature phase lag 

Jl, —-angular speed of body axes 

^ —mean free path; decay constant 

fC ^ cost/ ( 2^ rTav^ ) 
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TT^ f €) 

TT^ KVR‘'^+8K^^/(R^^r)TL^ 

TT+ “~“A/jf^ 

-’-'“‘minimum temperature peint from the twilight 
© — ^position of strip on the shell 
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LISTINGS OF SOURCE PROGRAMS 


The Gomputer eodes developed for the present 
studies were all prograraroed using the FORTRAN-iV 
iangua^eo Except for the Optimum Spin Rate program^ all 
codes have been suecessfully run under the WATFOR 
compilero However the present FORTRAN G and H 
coR^ tiers have shown difficulties to accept the codes 
which are usually not optimlzedo They contain some 
redundent statements left over during debug~peri od« The 
listing of the codes and some sample data are included 
In the original thesis on file in MIT Libraryo All 
codes were punched by the IBM 029 key punch and they 
were riRi on the IBM 360/65 processor,, 



